


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1999-06-01 


Examination of time-reversal acoustics in 
Shallow water and applications to underwater communica 


Abrantes, Antonio Adolfo Mendes 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/8359 
Copyright is reserved by the copyright owner 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist Ser Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ies) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 










































































ictal atta eral ite lil len liana die Nl lhe tthe hen eee litt te Be hd Ade inte A Bens tr fates shad a! bs HOT Ee eae am Sion 2 r oe he € 
erm ant eee otek to ee rec 2 BV Goad) wel e-gere er case + pegd pt ee eon, othe ge 4a Pea tawn tard! Saweoba iE 
“rs pe rwiys Tritt eur te amd A ereoa t Saas iS TERR TEC ge ented ELTGt ao enisenctesee Aicece ie 
cle inh okt elas OA tolerate pei Lee FE ache : es t-te miele 9 Poh ieee tie tye ep pbeberd lure ete te, 
? Oy ie Ps Py edt ead AoA etry A m= rh Bearer peace peer ot s Feepeeety Bt 

« t ids i) pi ey rire 
ee. Ra SEX ae ee eer ey Tipit ro) eee ety, eer eet renrar eer eres Terk Pee Pr eran pte 
are renee trae re A cio bball Spt te ttre “ splay pesto pan a . a} he Haalote st we iit tht 
et ae oe | SOG OO ie cela F wept es aie st Here Tor ara ee yer et Oey ag 





yew te 





agueqqes KR avebd 








Cerys 
eT ere ett ere nn bidet raleblale LL Attell by Seeded ot tay Lod beh Ld de he oo Dero ph iebids ho Me ali 












































Toe te eer toh we ty oa bi -tadhcheedhd 
Pht a 











ae TT a , 
iste ee hee ots tek 
at het ee Ce oy bg 
ee 


H = 4 y ’ oe te tor Pitre ttt De tt ae a ee it eer rae Uekn eu, eteta eles ie pitlatiealies 
af oe, Bee it pe ee a ieee f ho TL eae “¥ Epa a Thi hie edo oa eS . aes eet LTRS oi a aS Pra Ce eth avy eI Te TTT See beth 
nage Tether apr t ey . Bare phate ee tory ntti it Neer the a Ah Sie Lehi me =o Piha ea are aati AR aN eT ashe et bee TIC: baht hae hohe ad 1d Gay Pl dibcledebbs hdd toned gtniaende. 
pdb heh te ba by badi bey heh orb reap eptinek oe ree Leth ah Ts etibeh la tenon iielh nebo be cette ah Lite hor EDL Cee dat oe Fates Heal rae tb ds , se retard L 4 Peer Pe or ee Feary pelt eS 
ae eke Oe eel ap ie es bed batts hineh pe idrten ele tibe tadedete LL Lett eh tea the. od IP J Losin lind ip aad twee gecder ye” “ ih er Tit IM Tia ey ee ee Pee ee Col poemeey me bx Pte eh wihesde Leider teehee 
5 i CF i 


ok etiblbetth boliobias ch bbade bab ok ibis : 
FES acess oid tne tei peerpiepe Rts pd rete fo Me de hit hh tie a Pe ee ert thee ee 
De a eee et a Lot When last Dh he taka bie te. baal kek he Bek te Ge fad 

bald ioeal heede et Nerpel, Ld whee Shi alsin di obs dieatcaushlaheak sirtiacmes Saas 


Feeback ph dnearaete | 4G bhebegite- ste ie inedete e ttt sree itr ar eC eens oo 
sorted bedia teh rbd Mb latsadei deh Dan her o- $84 ta ca-Ghebe Deine thad What kamen eon alt ae a ee ee ee Sn Coes all 
ht ndandt yeeeaeenihnditr sweep lint tor Pp Nirholniivetnthic unde bode bl bert ih ck a Lae Lee tate or oe 
ronda me heversitie tenth rbipebeds leben mo ob Leib tt ee eee 
To th un ad bh Real of Seth Al thee hateeherbipt rid ns eta Ll nt ete ie ee ae tt es 
occa De ebadt-a: rh inecbada dated eb bagi aha ah selena Sebel odors oh a cee 
oe ee sod eel be 3 ot Bh taind. lor eadticheter eda deefie te batt ptide deutdenieal elated Gre atl IL eles tae Mr | en PORPECS Seam ieeat 5 acsahchel, Nala bth 
ete a ht de be ee ta ee oe Lt hn aed eee merrvrieotie i Lot heer et tar trte titan reer tart 
var nt eat, th arin berhehnoaabde’ eed ecb tere bind tater Lotektarethedela Lak tk et enn herbie eL te ede bear nn eta ihe bie EAA EO er CU e Ct Sih mm ta et tle ri Poth eee te 
a lncrtiviotebip Sit rnb pearhlor buat fe9rt-tep deem Mochethd eb pialid- ere hebelee Bo he Lh he ab hes Li etl hbk le LA TPL PE OR oe Te Pain ME tey tiny s Pee Pe Cee ears Co List) Spey bh alia ein tel 
Fe apm peat cpa sheets ar eres TT Cl ey eal wheeling bap rete ta id Ryeraiee 0" py ph Baws BF yf oe 0 ate wey ges inh tedden Fi ARPSh Oi esl bile oUt 

let @@wen i < <4 <a: w ere Bash Lee ho te et LT Det 08 ate tz ad toe td Pee a ey Dt) A dh el tes, ll r ae Ph 
ie Heche fic de si eS senby top eeei br tri ee ee pter eaten tan ae aps LF Pala Ach VEVSNE vopn a lL ete eerie iti welt t eee HPeere ere Crt | Pr 


veers Perea geo! what 9 Cra) 

bd hd eee neh oh eee thee Sh iy tard ot Mt ay wae re Te Sere ZS : ea! tp ie heb Ae Ba 
het pelt Madde el ete Fil 1 od eee ts aT Ey Co FeUbe tat Peg we mgh se r t 

rey Ler Oe etireri ttt tee) mei t 
CD ie Pd at eh th ee et ot Te 
Re ee ea! 
7 Preps 


tA ee ee Th Th ok Ue ory 
oT we me rhe eT er) 
PO Sto aes er ad) 


































. ca 
Te en aay Teer) 
3 adel BerPesend se has" age 





a ee le ee ed 
tT heen bids ated ale ba adel ited Bok ny it doen Od 
tel el tet oly te teenie mt itt te tte 
BOL) he bet lath hc thaeel 
er ed th oy al 
Trey Pt Sie te hia ee le) "ehOsDhe ier Fc ye) ee 
ako ahd SOT he eh elo le Db lie ache i Bat Del dl, tele sie he 
Ch eee Re hate eels bedel th ed ee eh hed ONT | 
bP CAdi on ares, 














r 
Perr tat ry Py 
SUSHOD Spy FR Ele DC aL oe ody byes Eb les eRe 
Aiea a eee rae ty 
4 alaat iar itis Lia) 





atRead ges ve stent 
BOS Cet ke heal billions 
Wage Cpegetqris tains ge oak 
Trott) 
ur ad 

























Bi bie y ry Cad *ees HHNED et 
Se h@ ehhh s DAC COE te tr tyre 

De de ho, CL athe ne LLL Bl) 
aL TP oe Ld eR Lael PP Pe Te res) 
" ee tel AN boot eh bl) 








OTIS ere ate Pit 1) 
kL Pe oe ee Pere Tis a ae 
Chad ee bee ORL Pete eee 
Pr AP" gheh p tenagpese? Perret se) 
% apna’ coteg' a WT tek bere 
ar TY 2 Boge g er ase Ne wang 
PO pa care Lace t tgs atoh gs dBer HPPsLee Meee 4 




















Cie eile Peo rer rit ory 

ct ee ae i ters? 

Pel eee horrors | 
Beet ae) 

























Saye 


OPPs EeeA ri gR eta esse gts 
Daal tl de hes del bet 


et te Aa ee bed a ds 


i) 
aq ate q Trea) tit eee t er 
ee PT al, mh rh lel ee ed ok De 



































Pe a ea ear em ee ee ee eee eae } Y pina Prrr Cervera ath ee ati kee cer ee ev er | 
ey eenhinaphaindterhiedabl Minha | Ssihathtletpt te be-thens etoeralieliba Lik ties tte te et ; rye en) Coen) ert ers ¥ Feed crs Lol baie tig tel ohn 2 Lal Tl easel oR Saha de td, alee ae i Rabel SAS ba 4 = Ld debate ede des 
meyer ae hpae PENCne Rey Sait f) iat Phebe deny te eet oA pat iad TNE me Ay pata 8 TS Ir Ler od eben ha papit pathy rortrtio eens ents ete Ot oe ee LAM tL ee as OE ed apr aeemr ne 2) weer ren oben ge chia! iy es ttle ot wp fe eee alta 
teh eet ee od ed et 





On Prere erent tral t02 
eer el 
rr. bli, ie ot Lil 
or Fn bene cd aed peep nth tel 3 
Pee te Pe th les ae 
ce Ft ae eek tel ee aati wate tte 
qhest sas gtr ins 3M a9 545 Gt 
Teer eer te tiles 
at Nebebcabe-doeheb-i dar Wa 
Pre oh es he toe] 





er i oe a Ae hd 
ce a! wwe Th td aed Don Ly te a eater oT 2% | 
ed ed bin it nde tahte tien el te Ll kt 
ee ee eee eS eee ee Pat or Me 
Pe ee Led Sid ae Yt ett het | 
De Le heh eel metry heir er 1 Ll 

aot} eet SPE PLY eae 


LL hn Si ee aie ere 
ie ete rs ed 
Chae Gia he tLe te ee ot 
SWE WETe nod bt paniy’ yale 
Yio ad ut chet ey Pee 
VAY Ter eee 
Ca oe oF i he 
Ca te to] 
ae 








PNT aie t tegen ta ee ure ys os ee ee rhe 

lel bade Rte he ee Soe it eC oie Cecks To et Peed 
Ce te) eh beth Aik ah Sodel PS Dead 
De th ee) kee Sted 
he et Dl Oe rt Cd bee Dae de 


Se tepgthne® a) PM ara toes OPA Det ee reg og th! 
Pe he) bebe dile edit det Le Seb teh det hy ad oor Dok 
Cr el eee ee Ree Le ere Poe eet ee Ted 

ET Re De i 








AL sat, eh beet b-dbadtin tod dk be hath at D 
ee te er eit) poe bet eel LY at od 

Ts De nh dteteietiealidih > aeaptiie Sealer ttl hatte tekcods bi Let Daa | eee 
Pad thelial) b Lelia tein Mehl tt elk Ct ree Or tre 
sa beat Lath tebe titemaietit ch etn alice tar tik Meade A-tehemd ae ha haba ceeds 5 SP WSTE we RAeeRs ati eaia- PR GPT LT 
cae heal he ti baepibaatr A-ha ital Ghee eae AD  theledh th tenet ol Alda tell et te te ee ee ty Lot eee ri Meee PY hm oe TY Oa Pee : 
eed he Tht ne The ih eT TL ee oki ei et tet ne Cet ts te ean a 
tr te SA ee te eat Ct te ere eee 
wo? aU y tle wore sie en abaftne Tes a OP ae ae Re dat eh 

St! 





























































eeetny Oey 
heir aay. | 
wt ahz aytar cette 








LL etalk bode ah td Se hel 











































Aig , owipe eet Se Para Lite el) trey en ted 
Re Bier eet ria Aer innate dc nis x Pra resets fs g dda epi ie toed Sh labepdon ett Stel ac beret bedi nbdeaely 
Palate hos hl ee tr iS ea een ete aeeee oe 5 OY ad oe aap SP Deb at lt te a i ee eet Pepe eet i het ‘ Loa Tat du dahh joa ahd AL: # aly RCA Eb SES sb gry [pdedn te Lb sane calieta aL acon daha hebenns 
at eee ere neat ie Nae rene oye ee Ph Pty ei iat 4 a . Ppa re wre i tre Proto eens one 5 J [etd tet ed Sl head a is sagt: re lmkp) gee ait wht pegs Arta eer ys Ps Sh 
a ie tot tes Pe ’ 2, Wie dc boc Beplitnah, Teton oll hn a alin on aitirin al labebep-tl vad) oe Biapeva casey PSs pee . Ph tat of i aes or ba | 
A et ces ee ee a eh I ak cal «44 H Oe et 
‘ ee et et at Se Seed aah tale | ered 


Pan’ a@ Oe gts eee 




















. Pret tetera! Peat Ty . He 
Pte te treet ait lato ts Ch i aa erred Fl al do 
PPM cy) eer hee Bed ie eee Sd et i ed ee a 


ba wages 
PeahD tal Peale thts 
oie er hee hk 





Oo a Se Be Se ee 
"A OTPRES E* ye, 









a ph prvi lfaemnllseeccheghepten Te he pied deethte der ee GLa et sleet eT 
atari chesah gin Mdecar k-le ea)-t dheptedete he ud teil kee ra er torr eer itt. er eee 





e 
bible Re LL 1 cen-e 









































f = i etry re eee) : 

- Letina. bn atic leindtanitinetiantiisiie etna ett an) kas salt ei eh od oa ie et OT Eide Rowe ghee Dk oh bas ed pa fat wre Pins oh er Lt) tu ee Lod 

cl ok cred ak a eer eaten. fe oe agar Oe Seer aed bh ke held Soo 8 ated lee kept EA) d cae Pier ie aT hr Ped 

1 de heptaaliatendd Ae phoma ste tehote oh hed aha ett ee Ge ee nr eee Pe : Seg Oe tt ee ee ee ay Te eee Layee eein: Sl ei bigest gle ala depos di Ita dhe ola bape Abed ee . ; La . 

Tah A leatlng ytd dedinets, wp tins bah hip ino ae telah Ot seis Me Lorn ee TTT thats eat ee Sy Bh bbe haleh chee Sh eld adalah atalino hat i talaatclan det Aa saab tie : 5 : ; , ar.%¢ tg ones POF eee hath penton chet eat 











sintphetharkapisterdliialh-tieeidh- deed th-tincia Iorda dhak cranked lat lehah Witte thot bee ee er 
Oa Le ibeatielis bh d heetaailss tee tal agian he ae Ee The a oak te dl 7 
cal teat tbe ertasitilalieeatshacetedetyltien delet yank eanlen dus tis -whs shat 
Seid oh nlp Ser habbit danetin beet e ee te al oa eee Toe tl end Pre 
wath detdaits inte edhe tanmet oe an tak 

ae el eh ee ee ee 
Stach dibey ga iactinnPnthathaattabp-efl-tonhetadeadenased lank Tie th badd hesitate ee Pitt ot Pr oe tt or tien 


ores Pete ee Cre att re - he Ge Sette ot peak oaeale pe Bets Py 
Saw ie Cte tia Ct eter eS trim eee Mee eG) J v FONE Te SORENTO aaa 
ee he tee te ti ee tT Cr iad Me Tier tr Cir kt eek) ey : TT on Lialhd ira sh 
7 Ce . Oo a , en o eer tn errs . 
Seat hod el betel Ls eed iol ae hes bh SL OA Eee ea oh Ly Ppt een ey Petes 
Pre et ene nos 
et Poy Cer ere re ie 
PA Re 


Ga PEE ay BHO eh w Lee Cur CFs Chaney FL ehiee aa Ae tt et tide oT hotel Tha he} 
he ee hh ee ke tata fete sane © 

an oe Ee ON ol 
Pada OS eg gY es Crea Oleh Lhe 2 
ee De A oe ee ee ee ete ee ees eP baie LED e! 
US bee ea ee Sr ee tr ns ey ee et eee il tO 
Sabebyiet res pete nates a 5 


ens'e fen 
oe eae CPP eee care Pre 
ret aesies ee Rees ae er elit 
Dd ate 
le Dee a 


























tslasalih i 1) bal hbase pierced tosh ek pd LL td te Pees ok on deel hs 






tira il ver est alee | i ty tt Cuan We ghtas 








































sf Pa us Countve eres] PPrr “ tt a he) Ce a ee ee eo es eee ee ieee ell ro et eer ae et 
sta Mat athena eek Mh damned ait bh reece eh ole Ll deus ut Bie Saat ee es PEs thio) per ir ecb taheh bP err plait bet pi a ee PRED Spinal yy shart Pa Dice Leta as eet ne andes DL lato on ak) pa ate Pes platted Fy ah a bo eal bed fp fo in 
hetip=jd tiered -shara-hentib take oh hie haan eam omnia EL ok oe eh ae eh the er ee ra eee Pyearetort 2 Oe lel ae oh la . Ls etoacth ea I te pened EL be tk dN teh et dal oe pashasy renews Te ed ee i ed ee ete el aed 
pinbadWhantirct-deuteapetiae Mode tot hang -Aocke lp tiecearadt 2 teath tikel nee hans oeh ted Peek Mabe let tere ee ed eee Se en kad DE doh Se : dine aa dR ch) Stabe 7 ie dh he Sten Foe) ae re tren ty ete] 

ied ocean sth Mant. 1 dh as a a hs ot Een te ee ie) Shenae! PL le ek eked Soros sana = be Ge wacee Rt, det te Bet i 


ot cE OLE CS OR Foe pe Shp eet iy May ar OF Ri yey Behe deat doh 854 Petedunt Feagd-e SERS ROR GSI” ener é erery irl Terr eee hy eee Naty on peOh Cat areas etee" yagi 






















FP eg reas 

5 ah - : Pr . 5 7 Bett aw LU ARR ASEE Ne r = TUTE ng ey = 2° hS2 Ty ir 

Shee itn eebagelr- Me Damper laberepnht ier teh Martel A lh aie Ue deck ok Un oi ok Be os ae ee een a det en Lek Le eed eat Pe hereatiiptl ch b-dae- dh Lit bel heii aide sih. tea Sel Tebe! pg hind Seca bleh ode ua Ct ah il a ae a he brad to) te. Sahel atleast Gi f eke rel LF et tl a a Pe ere Tt 

ee chad chetbcdek rate eet ee NEP Be Le i page ee Lette ie nla = H en ett Tree Tere tet cree ek es re er eee tant ST rer a bl eth tee tell Pe rp ped tats Osaene retutacs bY det Sel dal ol LI hen dal hd be alll Dh Dale el Tard & ee RT Reta eo tee BT Bree Pee er ey oe a er es 
pth cies Si Putlatictird doleaemmbablt wre ttr e tee ete on hehehe a Pele ree et apg et eae Tir rte | ST el tel bet ie Ri Rae CS 















B.deewg ween ef en a at 





: eos 3 S > Pa te or hone reg ee ee eT > Aveo ter, 
Cea te aT od oe ak bathed oh oleae tof oh ot 10 tal ot ide Tiel al 





PiaTees Ae Ban tedek se 











Cee eee ere. ance Ph itaet he) ot ee tera tl ir er tty ed Cee DD ey eel eel Pet ee eh Skeets, © Utes 
Te het dee Racin Pamichieed dit id aie et od ee et ee Rise apes ea haw: a loeede ene et tidied chp tscg’y dbottdinwldeabe a eile ei eek ar dy ee eee Ce ret ee ee Oe eh ieee he 
hat a taliteeteatinen inaketsiclibrattataalgte Web lath. ddhedicael deol a nl ta ee a Ck ie el a ce nee ee Pee et Cree td ote ee tt Spiele tallied tee he Te teiclnia etien sich Aistaan Titian, twat tahalaria' eda eielalted at tate Li b = 











oR ee eran ee Seber eet ee Sen ee 





: : : i OT ed ened a aS 
ee ee be De Ce et oe Le tee - een ery an TET ie ite art rt eee a EL Soe tele SS Ch aha a chalet Neti a bal dace Mt els i ell ade Si Da ada Bal bakit tla, Soles en tale 
ert say 7 Poe ts ie Ded ot Mel Lag al tin re Certs 





































































' 2 Crt] Add bh bl del ida be LS ad on te nr Te ee Pee had Sate ee tet ed ee Me ec en een 
Sitiah-thagh divah seats: ettheabetl da haiaibdiom teks Aath ae een ee ed ae Te ea oF oY es CE eat os tae hed atk Ltt Cr ts Oak nO na ee DO be ala dlalala ied chat ala iade 2 sete ey be a tt et) ie ee hee Lark hr weed SWS ap end tReet hy td tahe ts Ce eee ee ee er eh oT 
he tnnihd = e- “4 ry eertreTtr hevee Tree eer ever iP teres ee ont ors Tet ees CS Dahil Ah bed bal acted Alin an A deh thal hath i a del ail dad ache told Loot | bel bP aad Pir tatty Oa Mee LL ars tae Be bask tae en 
aad all A = Pe eee : Py ee Me ee ed a Phe hd eh bata ies je edhesthheh dette hi bod h deerdet Seabed ater id ee ee LT ae od P 
tat iietaeddinedl-thirtnratipainday 1 htt debates adikamtamied lene Ce eee eae ene Pe ets ee tat tt oe ta Sy Pe ail Lb aedbeh ar betlcedtia tal neat A! bol) Me Ti ecliies, ealt sat te FE OSD ge ZNPne "44 F082 
pndytntiaiideamdadittatieth bio ation el Pa ee pee fept Aetna ate NOs Oty thee nay toe ‘ rl inlet taba cies tn bel tas hadi il tn hahaa alia tle last. | 
Dieishetdeatiy statieieainalattin distin te tn dik betas ti ah ered eee tae tr) Spe ly epee eet ila seer ee, eee , be ah Cop hae ae folate det lta Collie : 
Se 6 OM el ans ie 8 Se an he Wh ge eG ew ee Te ee ae Se ht ee eo ney poy ee eet ie tee Pett Ty lel ctadeetiate dat taal Lh oe cakes eT SR ee ol Ell del delle alah hn dp gilda 2 Hy hd tig Ch at Ne ill deta 
than theta taht ape ch adhd he a Reet ret ot tt Seer Th a ie tka Pe ee ee eee eee ee Ge eT tee PN Cer etn bhp Mia aaah ore. BS dae ie 
fine tilled teeta tipple emasbetaltay* dea ata lak Te et ie tee teat Pitot ela chee to diiehabitadaad UL eat en Po Ld aki od (Coders ty Gas OR DO Om tly Re FP lg 
Seendde Renal ik teiiahattnentionah mudidoe dah tm he te et Pose ge eee ‘ Peete et tae eae pd elated asta iy emerald 
lbertaetheete Aietend aialiieal al inal al 2h ott cele he) lett ek ote) ie eae et eh eee eee et ae ~- oh a ted baad sta Taha fecha belated tate hd TTL a ¥ he a 
teeth teresa Sespesthdlin ha tty a caineath-eage-Adaetltaedi-tnhe al ite he tee hel ad TT eet fon ee oe Sp SE yy Oe ES pea “sue: ut pent of, TE tL tel eee | he Tak Saehdle si SIE ott MN bad ‘alec Din de. seahntin LLLY desdnahendh >. se hit Quai tatcinle tee skidded 
bathed ced alee ale ial tattle Manchin Ak ih eet Tied ot toate Cet ote ee et eT Pe WES yng ree ig a BeRU EHS A wena egeye = T"4 Pl viene Pret, ees ee ee ee Fe a Sal Sabet oe taal, alas ageeharhmn it Nah al ale aaa dia erat 6 maa laos, clan coreg baal eh haste ree ed he ae | 
htletlartlhcin sen addin leriaden beneath etieh tae ek oe he ae ee Be ee el ee ee eT oe at ie eT me tl itt Met eer ee wee or ee eee ea ee dk ne plettahdtt ted Keb deeb et Ae i ee od hl ae doe ad eer ere rh 
REC ES sche lth leper Set ee eee ek a er et eee en ee Ta Pye oe ety ees wa i = er eee wre eiee™ Arstve™ PyueGes*S) 8b wee aM OPER ey WE WOR ry Wile there plat pemk kd all eee AO Ter Te Rr a id Te 

SN Del ee eh ks eet tte ee tn ore Pe ee ae) ; : < eer An , ert ty Crt tt ee ree Tt cere re bet dpe Ucar’ ms ‘ St Medal Ds adel Peedsgueng 0 . Jeng Nae ag Sad y Sle ely om ter ote) et 
Oe ee a ee eo ae Pe “ ee eer ta Lee rere th ek ry ie Me oe od Re OEFe Sef FEM WE 1g = Pen eae. puasinns Baer ache Le A Sk bd oe tl hs) adalat he Ble die ete Sig ee ee ee 
sete teatldoalh  h nthnet Patiala Aat dE pieeate alt cde ie, De dee hh teh ahs me oT c PLP aT ey Tee et te Oe Cea Lhe ie Salih beieiniia ciel I ak inte een A Se Ord a ad ee he te ee te te 
hala ob batt ciate tte kel al Pade Eek tht ok ne noe ui pe Sas 4 [ar ps PS ee ee oe) Sh eee eT | 

LD ae 


Ca oer 
hulle alk pike eeslenditeriaet neta Teale Wt rate neath atedll. eae We ne oe a, nee 


Lacdsunabcittietduntionn’ qidediend meocumeks a nell atl ale ta hh et 
dep Un tiitibelyin ale, ahniettedonl aleeteatterlipeiel thle ehan hd tk ee te cn 


Core pe re Ce ee 


a 
ee eet ot ted ok bh) Le i, Pie od 
an addeen en 


ee Loe ee et ee) ee 
Stl Detelindigeh eta Biles Sead elie ait betel toes nn he 
tte i te te eke ee ee ee Penn ee 


5 
eet et Se eeites he hte RE Lalo 


Ce Et ieee ee ee oe oe i eC rT eee ety ere 


EE ae ee 























































* rr er erin! Is age a i . é Shoe ate bent drat ees Puneet ed ee bali lthes Wiel ictal Mid od ta dees daiae ales enlae. 1 aed a lieadandepoeiat ae ok se teal a 
-ghingihr atlanta sn cate theiindlle nthieatie aeal oh or te eae ee ee eer fom hd) ated alten hein inde: set de iahelial te he bet teal Sel Ad, dt Pamela J i he eee ok Poti i at ly we eee ee ee ee a oor ie 
oat edigthesthnatiehn de geet dr red Asa Oe eed Oe ee atk ee et ere ee Te ee Poser te each ggr ese he ote hah eh ee me Mant alti Lal SS Bette lode elas a aaa 4 s Ste tgs ‘cane vel oe edit Lh bd ee 
ee een  drti-iahatteh-te-Cude leah otal oe in ee ee Od eh fie ee bee o Gente Se eES™ = mee ngs dal Be A shah eldetheah aiid 4 bt et ne bt, deh bet de be bel Set ie al ent x eT te Ae ete el ea ne Pe eee te Te eT Pe ee Et 
etedieait ah ch) teeta, a het he erry Pi Meek lal Mccann, ae nd re a a ee ee ne Pe er ae pret te Sh I tele tel ated atelier sellin sit Wiech peli ke toed Del ed deed De Ls ee ee a ee ek eed blir oh pes ha, ef eee ee etry Les Py Ce ee een ten ea ee 
Se tiathe- tA tea tr tlebda ehe-e) ee RAs on ete Mie yg OTe eet? OER Si Gane pte 8 Mer atycenys hiao™ Beng OC a alate T ali od hd pinta iar ela dadaatintn a tin Sail hates tthe alist 4 bbe het Seditel dinde Aetna bbe bade tet tdead ih Led te Lette tate cei tel ee ae A Se le ee ee eh te he ey ee en 
Si cok endl eee tae ie oh i Sk a oe ak why ee OT eet ae eh Wet odd alt bel ee eh ee Pe oe ee a pe Oty . Le Sete Den tl etetiaiehtitiel de Sth te) ar Te ded the papenki se ieee ei Be Oe nie ee Le ee Set Gee te ey oe ree ede De oe een) SS 
 onnttieren ee tata atte athe te eel) bene LE he Le hae thd eh Z : . Pipi. hit ta ae tad tn eae 7 . do o rg? ad sete -lecceiaintie Libel | siete ne tte ah Pema ys “etn? eh Ae oie ote th dee ee eee SO te ee ee ee eee kt Ck ete Poet ee rrr 
Sentient asleep a ala ol a ed he late Ee oe ets ht er on tees eit ST a a i be ble th Cdn Stet de) Cee Let Or) LT Oma Weng Woe rit eh oe a Try tC ee eT ee Tes eres, 
vcilbp-chhaaiathh site-thedtiliehe-tatak-Aoniaekeel ok oie ot on ek ee ee ry pt En oo a Th PTE Errata ins ae Pee Dele he bel tetstel Leth i) a eerie ers Are te er er eee) 
- le ed * nt bs P z at . 3 <t pee Paar eT ee i er eae y “4 ee See an =. a — c o 
a Sia alt Sal piers tein ayiteatp Tk Daath het et al tel PR ede oo a ee . i Che Pee Perey P ore seer ner Pacey ee oa t TO PoP e OS 1ST E lel taal 2 ote PP ream Se el a he ie tel ie Del ba ellie! (ek St et ey Pte 





id 
a hal ah oe tet Dad Oke a ar | bat peer e 






ATT a taaWe f sea * ad rd Oo ede eb dle “ines she) tea Del bel Daeatel Midd del bl et ala ST SF ta & 





TT ie tr, ory nt 
sr ethentiahtneds, dhateette Loa alii at ko ki toate peat tated hah adele el he 


ie Se eat rg etre eee er ay 





ce 
Se ee oe ee eee eae Bee 



























bell peaeaaress 5 PY} siete doh taht rt ae] edie Sit eee Lie? Lie OTe on Pi ee a D 
meiepiincl por Neeti thle deel oliagll Lita in ah Mander cael tk tee ee ee rs ee eee ory Perky prec tlad tay te iy 2 a LOpt eemeaeg eX as bee eet erat? at) 
Leen tires tin ttanstinsehettantttaitn attend aie thim thea Tok hat tele eth ko ee ere Ce ee et ee ee ee eee | Repo tte See en | Ci Sh ST ta shan intent tached 3 nae Fe ee at te titties date te be 
beth ithaieatint- ie tiadatintiiath A dde antentein C ee e Pte ot ut ate et tt ea er nga yee = 2 = E pr ear ad ancsatqa iter as al cdantee cath eae tak hai aS “a "¢* COLT SE etl Dledal Dat eel et he he te PC a at eee yet 
siete ed caisadncali tie toh dp che drat Dodie dele See ah ot ee ae ee ee Cr ee er YP ely a ore eT a F ohentd anf One renete Sed elise wag te ie bie eee rt Cree eres tol - 
titi th atta chateata tet heel eine a oP teak ok ent td ee a a: = e Sn Suir - Rar = _ Perret re a 7 het eee Oe ene aT ty athe Se ee ee Se ee TL ll ae Pere eae or ere) Peer gt eee Ce fer i) 
ieclapaiieda dlaheads- theta eater, sethantes tietiee I ahaa et ch ieal th ke See ee eet ais sea ey eee eres et oe rd es ee he Pe bed de oh he ie en ee ie A ie aT ee te oot a Pi eel eet, Tie) Pe Py = ee Whe chia tay RNe whens 1S Og hay he, 
del hey lieag lich hetlienttin, tate eats alate ell ahead Le ate te on el ee ety re eee iia = aie paren ar eT rn er ah n° of Set ald bean Ol be ae | ee oe Yeh UY ee oat ee ee ed ee eee ee ey 
aidan te tt Wi de Ton a kod nebalttiae thal el a at et hh a do tk RO en ee ea r ta Pee en er he Pe et sk a) Sei telecine Web babii diet ait del Mael el e | A Lie Be “ DS ee ae ee ee ke eee ee 
byt hngntllrncindl-aptiatilienthe pa tatters lashed drakah i cdthatt th hell eh ae ee oe ee a . arr rrr. Cl ie at etl > “s ot MR ede 2 ed OS he teeter rk ody ee ee rr lie DR had dee 





hecpidlnntiiastinet tina tation bchtend-oth- site tirthas Latter tating 1. ahd ie ek ae oe wt 
a cthershrflleelnehnadinainciainatanlinn berithnatathn patient te beth dente at a 5 
Taphcken cttcthah—chag ainthrabatipaesinhepttninteaeah tnalitte tian a tee a ae 


<»* a 
Pee ee ee ee , PS a ne ht tbe fd A ed he) 





Ptererore metre tr. rrr tt Ld be ©" ee, A ee wet de, eh let he ee ek od oly 


Tr ee eT beter) 2 Cet Al fee te 





de hel et ke RS eal et ie te Eee 
Ch ee et ee eee So 





“mo ee er iy) ee oe ee | 



































































































































~ ey Me es To Pies ee eek ad Td ba et er eee eee ts re ee he ee a oe ee es 
si al tthe trata rhaalinch ht dieflaghneluelntaatlnatl-alet at Sanniod th ok sk eT na ee Yee Cre) a Te or be 5 Ce eat 3 hs beset ital Lal tl a teat CN Re ie ee ee Pe 
euihydp bap aetigelagteciphaetah aindathactedaseel atk adh eke ee fe payee ae re er res ee ee ‘ aL oNeden oer) eee ee ee 
ndiadiraitiedinadn teAiehetelaelnd ak-te tnd Mab th Tt kn ee a eee i ee eee A heel hen Ltd te : leptin te et thedd ora Oe ee ed ee ee en ee ee ee 
Saticliartetaaansincatiteled uth sh he Ad atl od ae Coe Swe asensto eo eee ea ea rae re ee 3 eee Cs dei lle ee ee ee 
s iientaatt- dike tanhaliarttiethe beaten Mt ak ak ot oe ae % CaM. od ol ad Agen am rare pins eee a Me cen eT ’ . Fi TOF FER iy gt hg a Or tate pee te gh ie te rol 
Sear Sealietanat di tis, tetitanieertht atid Did ate de a ok tek eid eer) c ¥ Ser tpn nee eas — c . PU ee re ee a Le ee ee eet Pe eer ee 
Sel pide th tales dt a a eek te ee ee bad eh ee eed a Pe eT | Se) iT at tel 5 hae Mc iad hte dei LL ae Se te ee ee ee ee eee es 
a Na eh chet deal the A tn ok de oa Fe ei aod Dd Se ees eee Te ee rer eet. heen et pu dhtidMah set tence L) OO FernPare yee ugh ya o ert Pde ae ee ee See Le ee ee eee, et ae eet oe er 
Peet eee 1 te SR ade LOE ene ae err ere ah La laledhed lisse) pp et i Nel beta ie teed bleh at eh beh Ld etehdeltialatel, Jb bh LOE. a > ceor : Ame ee CNEL ae ha eve De Ce eee ; 
eh thd an ne is 0 le le LOelh wt nas wlre ee " Pat . ¥ a 7 eae CIE Sb ae ol ete Le De bball eebtaliheh ah dell ok Diet Ss beet td Atle 1 ade § Da ry keg Se helen de at dh Ae et a det ee Le ho et ea 
eae Se ee aren 2s ‘ ee ee ee oes = Jt Pe TSS ar ah | 08%. Oh tee ia hee Pee ey te hie een ren Sed ee ee ee a. 
: eet Neal he neil f " we-= , oe are en TP ee re — ay Perr . Pe) 4 Ct | Thr i Na ad Son pV gate cvs mdmts~ fa et elite ils aie tee el es ee oe tO Te eee 7! See ne ae tt et ee Le 
panes natin act iteleal theta dll, chntieebetiedees dh nie el Rete Utah ah te i : a ‘ - a4 x ae Peete a er . CE tee Cee ts ae Bed LL ; he ey ert a ee eee Seg tre ee a ee ee eee ey oe 
eee bliin tinelhcaalie Macca te date the ele a tt od he es at co ee CH eee e = eee a gokdyes Pee ye ee ee o~ "ae teW, anh” oo, Pr Stele eben a ks ee ol ee he eh fe | ct | | anne 
wirciaplocliagind anlar peter) Meh es ee a OT Pe ee eee eee fe 3 SR A allel dite aldo Steed or 
ieninanion’y-litieatn dictate iin ttin sea tee _— ey en ee ee peat I arene ee ge a ape rey awry | a rr ee pe * i Cet) le tte Tey = 
Me aw oe & eee > b idl ehha dl ie te eee oh ee ee ee ae Be as cher eo ee eee ne te ae PE caer : ‘ PEER aT ct Pret ee oP "wha ‘p e~ CS Se oa) ee et LT tr ee ee rpenir es 
Cycle serps Ailes dol hoa heehee ale baler bhp hah tatch etinge tes yet ey See Ser en Nats aid 
path rdihelaalite hates Dect lit tek ea aon ek at Se he en eer Pa as aap ell Po ee ee rf a Me det Madi attal Metal tate ted ot he, Rates Bele Deed ete tee 
du ncnahimiean tall ciindeshesl Token a eee iden tian meenGin mauaele = « Pe dt oe . © bw eee tate eee hee tt ee eee er ee ee 
sat tn natty alain taken okeRecahe” Tht ee ee era aun ees TT ee gt meme eg 2 oe aida et Fc bidiindy kgdb re ee ee ee 
hae el cll ecdondibedlarlrth atte an putntattn dade adh, Bi tre ra ee ee ee nn a orate pee * RT wnat e= . » aoe Lane Cahdeteti “L ew, & ek eee | RN tin dsl ae betes een tendo. 
Sata titel rir cha aaa ala ie et le ie a Oe oe nd of thon ot od Ce ee a ere Prt aed y = ae a eS ee = eee a +e or ey S ial he oat te . ee z= Be tense len ty PP ee ty | eet ee eee ee teen ee oe et ond 
bdr tien tea tied htt Mthealedh Putas attest ent eee et eet eee i Lg seeun ett ctree tee Pre Sta F She RL I te sit Del hiatal mesedte g wear ot i Me bed les a! de oe PRE oe Oe te od ee ee) re pe ee ne ee 
liveth: Hemnlethaatle liegt pie mtn LO SO ae nae ee " PS SSRge ern os Ua ed eet | ae J edhecdhentia teal | le beate Se TRE a ir et te idetdie sede Let Le ee ee p te eae rene of ah Se: 
5 Ce ie ol CeCe Ce : me 0 ee oe) et Ce en iene = yattn <1 Ce beet bet ee ee ee eS ee ee Cae 
bees ee ee eed ok ee ee atJccmw" =» Wedlees Ye a de a A) tal hel Le ng rar Wee Pt Te = at oe ee ie a i 
Che ne ee ee ee ee erry rt - @8f » ge Se ee Sh hel all ot BD Deal tad ee ee ee Se ee ee pe ee ee 
ST ell 228 &® te ee ee Pt ee elt ary SU ke LE bathed I ad eek Pan” eens connoted aR Set “= ere oT eee re ok Toe 
i ee i ee Lot Raves @enneta to, a © & waged: 9 er hPOl ems ores arya es a ee ee ee Te ane) : Sere 
- em by sapien = Se © 42 Po - cel titintial th tedden ft ee ee Sea por Py Pay are ee a i of Te a i. eee ee | PP Oe - F a diated 
7 as iE aa > | cerehapapteg ppl ca! bliin ie LI hie Mh ahh Ce ee eT Le “~ « rn) Pee eer Sr re a . peop tet CTR Shel oe ela te Dea eek PEE pel oe hd ee el 
- bethenta tiie alata oh Resta ah a oe a at ae P a en ee ar * - 2 Pr =<t = = - ee, linet ra Betet ese tga vs ey any Rayeets! Se eee 2 ee bo ee iY 
aarti al oles la eh A, iit Mite kh he ne ee ee ot re a | - ea ee ee as alee 5 eh 7 @ weet s edd sd | gate Petteeedetdiaeda tin tl a dl alt Te etd he tn lat is tele ete 
lt oe lated el Rel ke eet eT es Cs ~n © K ae) Pe ee ee ee ee ee ee ee Te ee) ee « S a, oo a Se ee de ey tee he et ee tee LE eee tk PO oY 
Ce eee ee ee ee te et Sed dE Ya Se ee ee a). eae PA ~ eee 7 ame a. a ar Dh ee te - all Aaah ted a 2 lt Nee el tN Aletta tel To a tel ited leh ota et tee th he le ee Ped Ok Se 
el le elt ah ee eT ° Cee eee ee ae 7 e ee 7 ial Se | Se eee fe fe eke ee baila teh asain tetetien in tel tee teen aed 5 A deal lth eaten atm 
Ce ee ee ed ee ee Py av teks ae ee ae 7 Pte ek PEP Pe tl Lee Re ee Se bee ie ee ee a 
ek a ee ee en cP ot “Meee oF “= % Sar eae a Pay . e " % Cel tae tall Ce 








Oe te ee ee ee ee ee een nere | 


a he, ee oe a ewan Pe ee eee ae eee 


ee as 


Pe ee tee tt 
et eee ett Led : 
Pon tt eho aT fo Los 






© SWREe RH eee Fe + Me mary ® . 





ar é Pee eee 
Pees ea se aad 





ee Tien ad atin 









dae eo rt ye 














ey ae eed bees hh eee ek Po ee eo ae ee ere 
aa Poe. See es ya Ore ye a ae a nearing et ies “ =% « Fs , Pe ad iil. te Se cgmyp eB tt ognete ante & ee ee ee er Pe eet ne eee ee a ee en 
Lr z Se 5 4 5) md) a ar Sal ead Be ee Re ee et ee eee 
th ee =o = - Ce eee ofA" Cd Sree Ce ee ee ee ee ee ee ee) a ee ee 
eh Ea aay J ee Tal OO te a ee ee Le ae 2 tee ene 
En or een Te - 5 Aneta at . 5 a al et a he le ee ee ee ke eee 
=5@ er Ory - 









ke 





le aie Dh Re ee Re, hs dente ed te Ne i et ee Le Eh to ts te ttt ie ne 






















































ya oe Pere) « a ie eT he de er he oe] ee ee Le eel eed be kk aT ee ee ee ee ee 
at cM ol ee Pn ener n ee a, an ee ee Oe ee Pde a tM eee ee ee Se ee 
Oe ee ee ee ee ete er] Pe lite Mteetebe ca ee eee dy ee _* ee ee ee ee een Le 
Oe ee ee ee ee ee hints " - baled TT ed elie tli o ete ee eT he es eee ee ee eT 
-_—= Po on o ~ Pt ll alain 5 bd ol Ribak MEE Da tale Ad wre Ore re Oe ee ee eee ee ee ie 
‘= oe ed al ois ah aed salsa ees Dobbie eh st et eae oh ee ee et ee Le he oe te EY ee ee ee Lee 
et - wae “= « ie LL te eet oo. ee Ce ee ee ee ee) a ee 10 eng Qt tQ eee 
De ee eee ee % en oH 8 « ha Wal tet ee aa ee ee ee ee ee eS ee ee ee tee oe eee ee nen as wee 
Pd = D Seale) oy Le J ; ee lS een ee ee ee eee ee 
6 Sl eh et ea ee on ‘ eee bie I ale a > = tt tle Bi et eee 2 kr bedi thet latit Mn Miata deena Dad Th th) Pe ae 
Late a ae ho Per igus @ = lial landctha adil ee hot Be et ee et en eee Te) 
ed PO pad eee eeepc Cen SAAC aaa— a. - . oS oem sae ” “ - YY er ee) Ny ed eh oe . Te ee ee ee ee eee, ee eee te 
a ae A iad 4 Peet OD Sd aie ates ell ae hee a) Pee A ee ee es oe oe 
> ae ee ee ee ee ae Ce os er | ta Lae ee ee ee eet a er ee ee ed ek ee ee ee ee ee et 
at A hath. c = Cha Si LS a diate literttie lin dtl tl athe a dlehe, Aehee te Tick y olathe att cede 1 a on ete Ce et ee Le Le eee et ee 
ty i oe. eer ee ps : = 20 Woe ial a alata etalon tl Quien. 2 rie his » C Pe tent te be LO) hed tee ee ee ee ee ae ee ee es 
a ee Oe = es ' é A . s S " ahediied at ele ahah saad al de Sie ial eer - E o~ ctaty PrOb ne -RaSong ge PAG. ae te 8 ae ei ee ee ee ee 
— ee tet Pe cola halted ae Py he ee et eee thee ate” oS ee ee ee ee a eee 
ee te ed Ce oes ~ awa’ salt gta » J LS 2 ger te Ce he edit del ie eee ee es Ne ee en 
Biteiaeinnnndete te dhath the hd ee a oe J bbe 7” D pitas ie ee ee eT ary Pe ee Peer ee ae Tae) ee ee ee ee el 
LLJ eae oe we ow. Opa 2 = A ea ee Sa ee ee et G ee ee ee ee Le ¢ D 
e-> - * a eee ee es Pee eee ee es 
WY tee ee eh ter ei rele ee eh a ade ee 
WO € om Mad eile adh tal Rati slab Pek tok iat te De ee Ee ee ett eee ee he 
_ ~ lt Ee ial te ee el tke le Le Ee eli Del 
r-wein a ee ee Se eee Ee ee 2 eT Cee eS Le ee 
© -—. ear . errr ne oat rem te tO Nee y wi eno eee Oe 
lates . a eh. eee eS fk ee a. hee, 
s . =e - we ee ee ea ee ee ee ee + ae tite Asien hel 
ed aa ath ee eae i oe et eed Soe tt oe tee 
9 ™“e= & be 5 ae. ae kee 1. ee ee Lee ne fee eo TP Pe er 2 ee ae ee led . 
be i. eee ae ee oe ee Tee err ST ee Oe en ee ee ee ee 
St a ee ees y ee ee eee ed ee en a et a ert © 
es es eet Pann ee I ery ree 3 “ 2 dase wey" tape e = 
, Colibeateet af Pe Cen, ee ae “ ee er Se ee Pe a 
ar BS oT =: eT eee ed ee eee ae ee ° 
4. ca erin died ~ i Gt dy teat) Sie Loe Cag ro ae 5 7 7 . ee ee a 
, a ae) « rd Se ae « oe aes ee ere ‘ tet aye ee eee 
OE dialed ade = se a Oh tees .* 2 eee ee ar ee kl w os 
= i heli! . Stadia! (tes Cede ss ee C x 5 a ee ee ee 
a a a] - ee eid Sat oe he a a) Cee eae tel Met ad ees Oe eh ete -—* 
7 ye saps et ns Lr ad am a ed ee ee ed ec sae hs leadetinet te ae eed a ee ee ee ied ae | aoa LY la ee eT tl 
j a a) taeda eta Rt =e ee ey eee Set eee a er ee ee he, ee Eyes owe pemet == a 
p my nt a - haediies ey Chem ed Ps oe ee eS ee en he ee ee en Se a ie ee a ee 
. = a a ol Ad ad he A ed eg og at OS en ee ee ee i re a) ce ee eh ok teed ol 
Ps 2 - s S . «Ba, ne hte te 6 et te he ee ed ee, eo sedeadiee deen on Ee Lt a 
i - ry ieteed ed a "a ee ih ee ce ote tot ee ek ee Soe tn Se ae 
= - fares CO ok ae ed cs . aad a 880 hae ng Mee ~ ee ee ee oe eer es Se ee Pe wat Ge oe ie 
ahs = A e P 5 = , = es P ear 4 eae S ae ee ry - ® rye ore a Te di et ee ee ee ete i. ae. et; ae Paes ee a 
= = a Bur be hat be a - . . a F * : > ae we 7 a on ere we Wee ee ots pit ue Date ie eee Cli thle Nl. Phat dle ote ae de 
~hacallied ba = a a - Tok es ee 7 . “ o a " S ee we ee - : nd Limbelhinbidinabediie Tinta tn pales A ba bade en eee hn th ol doealichetitel aie teed a is 
- a. ee a Ve ae eit I ase — bd aad li « eae a ad aN tl a data en all a sdk iether dled nalin atl aes epcethieadinc teed De Riedie, ainda Tdent - 
EFS - a pee Ry - E bad -_ 71. ad - ke eh ey -ea% g. a ed ee Ce ee Ee eh . ee) oe a od eh ee re = 
PSE = Ci en ~~. . * id enemy r , Lee Th od tend, oan inh la dioah Ant das deine PT te ae ital Rated te Lie a TD aw a 
Say Pepa nepal na a a bbe to i ES ebidedtdiieetat-o-dt ee rer eer nn Or ee ae Oe ee te 
< “o Se - = a@brasus os ay =e o <a 2 -« Sh md 4 ye owe a hed — eth) ad 
re a S u ro 3 5 * a = - - 


DUD...» aNOX LIBR 
sig ae SE Ox UE T" SCHOOL 
MONTEREY CA 93943-0101 








NAVAL POSTGRADUATE SCHOOL 
Monterey, California 





THESIS 


EXAMINATION OF TIME-REVERSAL ACOUSTICS IN | 
SHALLOW WATER AND APPLICATIONS TO 
UNDERWATER COMMUNICATIONS 


by 


Antonio Adolfo Mendes Abrantes 


June 1999 


Mines ce Nt core: Kevin B. Smith 
| Andrés Larraza 
Monique P. Fargues 





Approved for public release; distribution is unlimited. 





Public reporting burden for this collection of formation is estumated to average | hour per response, uxcluding the ume reviewing structions, searching existing data sources . gathenng and maimliauung the data 
needed. and completing and reviewing the collection of information. Send comments regarding this burden estumate or any other aspect of this collection of information, including suggestions for reducing this burden 
is We ashuvig ro Headquas ters Sty rmmKees, Dutviora® ter laonnaon Uperatois aid Keports, p25 Jessel! Las Highiw ay. Sunt toe, Au agit, VA 22202-4502, am iv ie UilKke U Managemen ang Budget, 
Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE (Leave Blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 
June 1999 Master’s Thesis 


4. TITLE AND SUBTITLE 5. FUNDING NUMBERS 




















Examination of Time-Reversal Acoustics in Shallow Water and 
Applications to Underwater Communications 






6. AUTHOR(S) 


Anténio Adolfo Mendes Abrantes 





7. PERFORMING ORGANIZATION NAME(S) AND ADRESS(ES) | 8. PERFORMING ORGANIZATION 
REPORT NUMBER 


| Naval Postgraduate School 
Monterey, CA 93943-5000 


9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING / MONITORING 
AGENCY REPORT NUMBER 


II. SUPPLEMENTARY NOTES 
The views expressed in this thesis are those of the author and do not reflect the official policy or position 


of the Department of Defense or the United States Government. 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 
Approved for public release; distribution is unlimited. 


13. ABSTRACT (Maximum 200 words) 
Time-varying multipath propagation is considered the most important difficulty in shallow 


underwater acoustic (UWA) communications. To compensate for the time variability of the channel, the 
receiver must use an adaptive algorithm for adjusting its parameters. At high symbol rates, intersymbol 
interference caused by multipath propagation requires large adaptive filters, increasing the computational 
complexity at the receiver end. 

This thesis presents a time-reversal acoustic technique (implemented with a phase-conjugated array 
or PCA) that generates a spatio-temporal focus of acoustic energy at the receiver location which reduces 
distortions introduced by channel propagation (including multipath), allowing the use of low-complexity 
receivers. Numerical analysis shows that for different PCA geometries (element spacing and aperture 
sizes), the PCA focus footprint does not appear to significantly change its dimensions. Furthermore, the 
aperture size plays a more significant role than the number of array elements in a PC array design. 

Specific examples of novel UWA communication systems utilizing time-reversal focusing are 
introduced. Current simulation results suggest the potential for high data transfer rates compared to 


existing noncoherent UWA communication systems. 
14. SUBJECT TERMS 15. NUMBER OF PAGES 


Underwater Acoustic Communications, Command, Control and 94 


Communications 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION — | 20. LIMITATION OF ABSTRACT 
OF REPORT OF THIS PAGE OF ABSTRACT 


Unclassified Unclassified Unclassified | Unlimited 


NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 











Approved for public release; distribution is unlimited. 


EXAMINATION OF TIME-REVERSAL ACOUSTICS IN SHALLOW WATER 
AND APPLICATIONS TO UNDERWATER COMMUNICATIONS 


Antoénio Adolfo Mendes Abrantes 
Lieutenant, Portuguese Navy 
B. S., Escola Naval, 1993 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE 
IN 
ENGINEERING ACOUSTICS AND ELECTRICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 1999 





ABSTRACT 


Time-varying multipath propagation is considered the most important difficulty in 
shallow underwater acoustic (UWA) communications. To compensate for the time vari- 
ability of the channel, the receiver must use an adaptive algorithm for adjusting its param- 
eters. At high symbol rates, intersymbol interference caused by multipath propagation 
requires large adaptive filters, increasing the computational complexity at the receiver 
end. 

This thesis presents a time-reversal acoustic technique (implemented with a phase- 


conjugated array of PCA) that generates a spatio-temporal focus of acoustic energy at the 


receiver location which reduces distortions introduced by channel propagation (including 


multipath), allowing the use of low-complexity receivers. Numerical analysis shows that 
for different PCA geometries (element spacing and aperture sizes), the PCA focus foot- 
print does not appear to significantly change its dimensions. Futhermore, the aperture size 
plays a more significant role than the number of array elements in a PC array design. 
Specific examples of novel UWA communication systems utilizing time-reversal 
focusing are introduced. Current simulation results suggest the potential for high data 


transfer rates compared to existing noncoherent UWA communication systems. 
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I. INTRODUCTION 


Traditionally, applications for underwater acoustic (UWA) communication systems 
were almost exclusively military. In the last 10 years there has been a growing need of 
UWA communication systems for commercial applications. As a consequence, there was 
a tremendous increase in research and development of UWA communication systems. 
Applications that have received much attention lately are secure military communications, 
pollution monitoring, remote control in off-shore oil industry, video telemetry (using 
camera images or side scan sonar images from remotely operated vehicles) and collection 
of scientific data recorded at benthic stations without the need for retrieving the 
instruments. Many of the applications being developed are now calling for near real-time 
communication with submarines and remotely operated vehicles. As ie UWA 
communications channel has limited bandwidth available, the bandwidth efficiency 
becomes an important issue fora UWA communication system. 

The shallow water acoustic communication channel is characterized by strong signal 
degradation caused by multipath propagation and high spatial and temporal variability of 
the channel conditions. In any underwater acoustic environment there is limited 
bandwidth available due to transmission loss which increases with both frequency and 
range. This is a major constraint in underwater communication systems design. The most 
important of the difficulties encountered in shallow water acoustics digital 
communications is considered to be the time-varying multipath propagation (Catipovic, 
1990, Baggeroer, 1984 and Coates et al., 1993). This leads to a requirement for powerful 


and reliable receiver algorithms for signal processing in a shallow water environment. In 


I 


recent years there was a large effort made in the design of techniques and algorithms to 
overcome this problem. 

In a digital communication system, multipath propagation causes intersymbol 
interference (ISI). As an example, in a medium-range, 6km long, shallow water (200m) 
channel, a typical value for multipath delay spread is 30ms. If the communication system 
is signaling at | kilosymbol per second, the ISI will extend over 30 symbols, making it 
difficult for a receiver to recover the original message. 

In shallow waters, multipath is eae due to reflections at the surface and bottom 
of the channel. The multipath structure depends on the channel geometry, frequency of 
transmitted signals and on the source and receiver locations. On a small time scale, the 
most important contribution to channel variability is surface scattering due to waves 
producing time-varying multipaths. On a larger time scale, channel variation is due to a 
variety of effects including internal waves, temperature gradients, and currents in the 
sound speed profile. 

The ISI and strong phase fluctuations caused by multipath propagation and temporal 
channel variability have in the past led to system designs based exclusively on 
noncoherent detection methods. These systems have poor spectral efficiency and low 
signaling rates and are still in use when robustness is the principal requirement in system 
design. To overcome the problem of ISI these systems employ signal design with guard 
times which are inserted between consecutive pulses to ensure that all the reverberation 
will vanish before each subsequent pulse is to be received (Stojanovic, 1996). The 


insertion of guard times obviously results in a reduction of the available data throughput. 


In recent years, the feasibility of bandwidth-efficient phase-coherent modulation 
techniques for UWA communications has been proven (Stojanovic et al., 1993, Stojanovic 
et al., 1994, and Stojanovic et al., 1995). Several systems have been proposed and 
implemented using differential and purely coherent detection methods resulting in 
increased data throughput. In order to accommodate higher signaling rates, these systems 
must allow for ISI in the received signal and be able to track any phase variations due to 
changing channel conditions. In order to compensate for the ISI these systems employ 
either some form of array processing, for exploitation of spatial diversity, or equalization 
or a combination of both. 

Array processing for multipath suppression has been used both at the transmitter end 
and at the receiver end. Coates (1993) and Galvin and Coates (1994), describe an 
approach that uses transmitter arrays to excite only a single path of propagation. Long 
arrays and careful positioning are required to ensure complete absence of multipath. In 
general it was found that this technique is more effective at shorter ranges. 

Another approach (Howe et al., 1994, Tarbit et al., 1994, and Henderson et al., 1994) 
is based on adaptive beamforming at the receiver end. It uses a least mean squares (LMS) 
type of algorithm to adaptively steer nulls in the direction of a surface reflected wave. It 
was found to decrease in performance as the range increases in relation to depth (Tarbit et 
al., 1994). In order to complement the performance of the beamformer, an equalizer was 
proposed by Howe et al. (1994). The equalizer is of a decision-feedback type and operates 
under an LMS algonthm whose low computational requirements allow real-time 


adaptation at the symbol rate. 


A different approach based on purely phase-coherent detection methods is described 
by Stojanovic et al. (1993, 1994 and 1995). The signal processing methods are based on 
joint synchronization and equalization to counter the effect of phase variations and JSI. It 
incorporates spatial signal processing based on combining diversity and a fractionally- 
spaced decision-feedback equalization with a recursive least squares (RLS) algorithm. 

The channel distortions encountered in the underwater acoustic channel require 
complex receiver structures resulting in high computational requirements which may 
exceed the speeds of available hardware. In order to compensate for the time variability of 
the channel, the receiver must use an adaptive algorithm for adjusting its parameters based 
on a beamforming, equalization or combined approach. At high symbol rates, the long ISI 
requires large adaptive filters therefore increasing the computational complexity. 

The aim | of the present thesis is to provide an original and elegant solution to 
compensate for the distortions of the shallow water channel. This solution is not based on 
sophisticated processing algorithms but is equivalent to matched field processing with the 
filter matched to the impulse response of the ocean. The low computational load required 
in this technique is due to the fact that it uses the ocean itself as the matched filter for the 
acoustic propagation between source and receiver. By using rather simple signal 
processing at the transmitter (suitable for real-time implementation) the multipath 
structure at the receiver end is virtually absent, allowing a reduced-complexity receiver 
structure. 

This technique uses time-reversal acoustic (TRA) arrays to generate a spatio-temporal 
focus of acoustic energy at the receiver location eliminating distortions introduced by 


channel propagation (Abrantes et al., 1999). Throughout this thesis the Monterey-Miami 


4 


Parabolic Equation (MMPE) acoustic propagation mode! (Smith and Tappert, 1994) is 
used for time-reversal acoustics numerical modeling. Figures 1.1, 1.2 and 1.3 present 
numerical modeling results of an acoustic pulse propagation from the source to the time- 
reversal array (forward propagation) and the propagation from the time-reversal array to 


the source (backward propagation). Figure 1.1 illustrates single pulse propagation from a 
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Figure 1.1. Single Pulse Propagation in a Shallow Water Channel. 
source located 150m deep at the origin to a vertical array 6km away. Figure 1.2 illustrates 
the transmitted signal features and the multipath arrival structure of the recorded signal by 
the 150m deep array element. As seen in Figure |.2 the transmitted pulse envelope has a 
-3dB tempora! width of 14.1 ms. It has a carrier frequency of 800Hz and its spectrum has a 
-3dB bandwidth of 36.6Hz. The transmitted pulse spectrum is given by the normalized 


coefficients of a Hanning window spanning !00Hz about the carrier frequency. The 


complete description of this range-independent channel can be found in Section 4.A, p. 


34. 
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Figure 1.2. Transmitted Signal Spectrum (First Panel), Transmitted Signal Envelope 
(Second Panel), Transmitted Signal in Time Domain (Third Panel) and Received 
Signal on a Single Array Element at Range 6km from the Source and Depth 150m 
(Lower Panel). 


In a time-reversal acoustic system, the distorted wavefield is recorded by the array, 
then transmitted in a time-reversed fashion (last signal recorded becomes the first 


transmitted). The time-reversed wavefield generated by the array propagates through the 
6 


medium and optimally focuses at the original source location. Note that the time-reversal 
process at the array is equivalent to a phase conjugation in the frequency domain. For this 


reason, the TRA array may also be refered to as the phase-conjugation array 


(PCA). 
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Figure 1.3. Time-Reversed Propagation from the PCA. 
Figure 1.3 represents four snap-shots of the wavefield propagation from the PCA to the 
source. The spatial focusing at the source location is clearly observable in the lower left 
panel. 
This technique can be considered a self-adaptive process that automatically 
compensates for any distortions due to propagation in the channel as well as any PCA 
imperfections. Clearly, due to the time-reversal transformation at the PCA, the received 


signal has the time-reversal signature of the signal which was transmitted originally. For 


example, a linear up-sweep frequency modulated pulse will appear as a down-sweep at the 
focus. 

This phase-conjugation process takes advantage of spatial reciprocity which is a 
property of wave propagation in a medium with stable refractive index. An acoustic 
phenomenon Is said to satisfy spatial reciprocity if interchanging the positions of source 
and receiver does not alter the resulting wavefield at the receiver (Jensen et al., 1994). In 
general, reciprocity requires a static environment. TRA focusing requires spatial 
reciprocity in order to construct the exact time-reversed wavefield that will focus at the 
original source location. In this way, all the PCA elements transmit back to the source 1n a 
reciprocal fashion, and the time reversed wavefield will have the multipath structure 
undone at the source location. 

Since the PCA only “learns” the medium structure that is imprinted in the incoming 
signal at the instant of its reception, it cannot compensate for refractive fluctuations that 
occur after reception of the signal. Also, if the environment or its boundaries change (due 
to surface waves, internal waves, temperature gradients, receiver or PCA motion, etc.), 
reciprocity is violated and the focusing properties of the PCA are degraded. To 
accommodate for temporal changes in the channel, the PCA must update the transfer 
function of the environment with some periodicity that depends on how fast 
environmental changes affect the focusing properties. 

The reminder of this thesis consists of five chapters. Chapter II presents an overview 
of time revcrsal acoustics theory. In particular it analyses the cases of time reversal 
acoustics in a range-independent channel and a more general static range-dependent 


channel. In these ideal situations a closed-form solution for the field at the focus can be 


evaluated, and some focusing properties can be inferred and may be extrapolated for more 
general environments. Chapter III presents an overview of the Monterey-Miami Parabolic 
Equation (MMPE) acoustic propagation model (Smith, 1996). This is the acoustic 
propagation model used throughout this thesis for time reversal acoustics numerical 
modeling. In Chapter IV the focusing properties of the PCA are studied through MMPE 
numerical simulations. It shows how the PCA focusing properties change when the array 
operates at different carrier frequencies, and when the array length or the array element 
spacing is altered. It also presents the frequency characterization of the channel and gives 
examples of how range shifting of the focus location can be attained. Chapter V describes 
applications of PC arrays to UWA communications and Chapter VI presents a summary of 


conclusions and identifies some aspects of PCA modeling requiring further research. 





Il. TIME REVERSAL ACOUSTICS THEORY 


This chapter presents an overview of time reversal acoustics (TRA) theory and its 
close relationship to matched field processing techniques. In particular it shows a closed 
form solution of the focus field in a simple range-independent environment. It also makes 
a generalization for a static range-dependent environment where an approximation of the 
field at the focus location can be evaluated. From this approximation the focusing 
properties of PC arrays can be evaluated. 

A. OVERVIEW OF TIME REVERSAL ACOUSTICS THEORY 
For simplicity, consider a waveguide with pressure-release surface and rigid bottom as 


depicted in Figure 2.1 where the axis from o’ to o will be the reference axis for 


Ocean surface 


R 





Figure 2.1. Waveguide Geometry. 


propagation from the receiver array to the source. From the axis location in Figure 2.1 the 


I] 


coordinates are related by z’=z and r’=R-r. In this environment the frequency-domain 
wave equation, or Helmholtz equation for a range-independent waveguide can be wnitten 


as (Jensen et al., 1994) 


[V> + k2(z)]G(r, z|z, @) = -8(r)8(z—-z,), (2.1) 





where k2(z) = is the acoustic wavenumber for a waveguide with sound speed 


(2 
c*(Z) 
c(z), @ denotes the angular frequency and G(r, z|z,,@) is the frequency dependent 
Green’s function in cylindrical coordinates at location (r, z) due to a point source located 
at range r = O and depth z = z,. Note that the depth z 1s taken positive downward. 

From the point of view of linear systems theory, at location (7,z) the Green’s function 
G(r, z|Z,, ®) represents the ocean impulse response due to a point source located at range 
r = O and depth z = z,. It can be shown (Jensen et al., 1994) that G(7, z|z,, ) satisfies 


the principle of reciprocity, 


>» as 7 >» » |; 
P(7s)G(r|r5 ©) = P(r)G(rs|7j, ©), (2.2) 
f ‘ ‘: = Ps X : 
that 1s, the acoustic field measured at 7; due to a point source located at rs is equal to the 
m : » , : x 
acoustic field that would be measured at r; if the point source where located at r; scaled 


by the ratio of densities at r, and rj. Since the density is nearly constant throughout the 
ocean, this leads to reciprocity of the acoustic field itself. 

The reciprocity theorem plays a major role in time reversal acoustics because it 
implies that the field from an array element propagates to the source location in a 


reciprocal way as the field from the source to the array element. The acoustic field 


[2 


— 


recorded by a given array element has imprinted on it all the channel effects including 
multipath. Due to multipath propagation the signal recorded by a given PCA element will 


have multiple arrivals. Figure 2.2 illustrates the multipath arrival at a given PCA element 


Reflected path (2"¢ arrival) 


_ . Ocean surface 


Source Direct path (1° arrival) 


Figure 2.2. Multipath Arrival at Array Element. 





where the direct path will arrive first and the reflected path will arrive second. After 
reception and recording, the received signal is time-reversed and transmitted back to the 
source. The back propagation from the array element to the source 1s illustrated in Figure 
2.3. The travel times along the reciprocal paths will be the same, and so the multiple 
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. . Ocean. surface 
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Figure 2.3. Time Reversal Propagation. 


arrivals will arrive simultaneously back at the source location. Thus, at the source 


location the multipath structure will be undone. 


Ie 


At this point it should be mentioned that the term “time reversal” is used in reference 
to the time-reversal of the recorded signals performed at the PCA. It does not imply time- 
reversal of the waveguide effects because there is no way to reverse propagation losses in 
the channel transmission. Also, “phase conjugation” is synonymous with time-reversal at 


the array because of the time-reversal property of the Fourier transform. For a real signal 
one can write F[s(t)] = S(w), where F[s(t)] denotes the Fourier transform of s(t). 


Thus, the Fourier transform of the time-reversed signal s(-t) is then given by 


F[s(-t)] = S(-—w) = S*(@). Therefore the time-reversal transformation is equivalent to 
complex conjugation in the frequency domain. 

Due to the range-independent nature of the waveguide, a separation of variables 
technique may be used to solve the wave equation, Eq. (2.1). The solutions are then 
defined in terms of depth-dependent eigenfunctions, or normal modes, and range- 
dependent Hankel functions. In this situation, the azimuthally symmetric solutions have a 


far-field approximation given by (Jensen et al., 1994) 


iK_r 
e m 


IKn 





G(r, zz» @) = (2.3) 


. nth? oo 
—i_{) Wn (Z5) En(Z) 
Pc. V8 2 


where r is the horizontal distance from the source, G(r, z|z,, @) is the pressure field at 
location (7,z) due to a point source at z, scaled such that the pressure is unity at r=/m, and 


Y,,(z) and K,, are the normal mode eigenfunctions and modal wavenumbers, 


respectively. 
The eigenfunctions and modal wavenumbers are obtained by solving the eigenvalue 


problem (Jensen et al., 1994) 
14 
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eae eel ye Kk ak | — Oe oa 
0.) a ZY nl ) [kK (z) — K-m]¥,,(2) (2.4) 


The mode functions V,,(z) are orthogonal and normalized according to 


oO 


ea) 
ae Oe 2 
| P(z) aa 
0 
and form a complete set (Kuperman et al., 1998), 
— YP (z)¥_(z 
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When the source transmits a pulse s(t), upon reception the a element of the PCA 


records the time-domain signal 


P(R, zt) = fav. Z;|Z ys w)S(w)e" "da, ey 


where S$(@) is the Fourier transform of s(t) and R is the distance from the source to the 
PCA. Due to G(R, z;|z,, @), the ocean impulse response, this equation incorporates all 
the waveguide effects including time elongation due to multipath propagation. For 
causality considerations, assume that the received signal p(R, z;,t) lies in the interval 
O<t<t, where T is large enough to include all multipath arrivals. In order to maintain 
causality the time reversed signal used to excite the fe element of the PCA is 
D(R, a T—t) where Tis such that 7 >2t. Under these considerations the time-reversed 


signal becomes 


p(R,Z,T-t) = [ co, z;|z, @)S(w)e OP dw. (2.8) 


Using the conjugate symmetry of the transmitted pulse, 1.e., S(-w) = S$*(@), and that of 
Green's function, 1.¢., G(r, z|z,,-W) = G*(7, z[Z,,@), and reversing the sign of the 


integration variable w, Eq. (2.8) comes 


POR, zp T-1) = | G(r, z|zp-@)S(-a)e'*7"(—dw) 2.9) 


co 


= | Gere | ee w)S*(w)e Je "dw 


=o 


where the quantity in brackets is the Fourier transform of the signal transmitted by the j”" 
element of the PCA. This corresponds to the signal received by the element after time- 
reversal and delay. From Eq. (2.9) it becomes clear that the time reversal transformation 1s 
equivalent to complex conjugation in the frequency-domain plus a phase factor to observe 
Causality. 


Upon transmission back towards the source, the time reversed acoustic field 


Prp,(r’, 2|Z;, @) generated by the ia element of the PCA is given by 


Prra(’’, Ze @)) = Ce Z|Zj. W)G*(R, ee w)S*(@)e”’ , (2.10) 


, | . , th 
where G(r’, z|z;, @) is the field at location (r’,z) due to a point source located at the i 


element of the PCA. The total field at (r’,z) is given by the superposition of the field 
produced by each individual array element, 
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se 
Prra(r’, 2|0) = SG", ei) GORE ZZ w)S*(m)er". Qn) 


a 
By Fourier synthesis, the ttme-domain signal becomes 


J 
Prra(’, zt) = ¥[G0, alg. W)G*(R, z;|Z. @)S*(w)e OT Ol ey | 2) 


=] 


At the source location r’=R, this expression becomes 
Prra(R, Zt) = ¥ [ar Z,|Zj, W)G*(R, Zs Z,, W)S*(w)e” e dw. Cr) 
jel 
Using the reciprocity principle, Eq. (2.2), with constant density p(R, z,) = P(R, z;), leads 


to 


GCR 2.|2, 0) = GUC za. a). (2.14) 


Therefore the time-domain equivalent of Eq.(2.13) is given by (Kuperman et al., 1998) 


Prral&, Zt) = (2.15) 





mi). [fs +r" BUI Belz, ,0)2/(R, Z| Zjs coat’ |s(u”— 1+ T) Be 


(2n)° 


where 8,74. (R, Z,|2;), 87(R, Z5|z;) and s(t" —1+ T) are the time-domain representations 


of the Green’s function and the probe signal. 


The following observations can be made about Eq. (2.15) (Kuperman et al., 1998): 


¢ The auto-correlation of the Green’s function corresponds to a matched filtering 
operation with the filter matched to the impulse response of the propagation from 


th . . 
the source to the 7 element of the array. This operation reduces the time elonga- 
tion due to multipath propagation, producing a form of temporal focusing. 


¢ It can be shown that the sum over the array elements is a form of spatial matched 
filtering analogous to that employed in the Bartlett matched field processor (Bag- 
geroer et al., 1993). 


¢ The sum over the array elements further improves temporal focusing as the side- 
lobes of the matched filters for each array element tend to average to zero whichis 
also analogous to broadband matched-field processing results (Bnenzo and 
Hodgkiss, 1993). 


¢ The integral over t’’ is a convolution of each matched-filtered channel impulse 
response with the time-reversed and delayed transmitted pulse. As a consequence 
this pulse is not matched-filtered because a time-reversed version of it 1s received. 


B. TIME REVERSAL IN A RANGE-INDEPENDENT UNDERWATER 
ACOUSTIC CHANNEL 


In order to evaluate the vertical focusing properties of the PCA, consider a range- 
independent channel where a closed-form solution for the field at the source can be 
evaluated. In a range-independent UWA channel the TRA field at location (r’,z) can be 


found by substituting Eq. (2.3) into Eq. (2.12), such that 





, S*(w)ee" Bini (2 Jet mci) hace een cae i(K,,1 — K,R) 
Prea(?’,z0) = YS y e ONS 
oo 9. 8nJ/Rr’ 4 PAK nk, 
j mon 


For an array that adequately samples all the modes, the summation over the PCA elements 


can be approximated by an integral 
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Using the orthogonality property of the modes, Eq. (2.17) becomes 
is @ iK (r’-R) 
Praa(?’s 2.) = Sewer — mer alee __ (2.18) 
re SNK, VRr’ 
Evaluating Eq. (2.18) at the source range, r’ = R, yields 
, Beene) 
PIR) = (oles yy ee 2.19 
raa(R, & ) = S*(w)eO" DORK (2.19) 


m 
Finally, assuming that K,, is nearly constant over the contributing modes and using the 


modal closure relation Eq. (2.6), 


Eee (O any i en ae oma eae 
Prralk, a) W)= STRK 0, eae 8nRK 0(Z-Z,)- 


m 

This equation clearly shows the vertical focusing performed by the PCA at the range 
corresponding to the source location. It is worth mentioning that this vertical focusing 1s 
due to the closure property of the modes. If the array has only a few elements and cannot 
properly sample the propagating modes, one can conclude that this focusing will be 
degraded. As will be shown in this thesis, numerical simulations indicate that even with a 
small array aperture, a good quality of focus can be achieved. 


C. TIME REVERSAL ACOUSTICS IN A STATIC RANGE-DEPENDENT 
UNDERWATER ACOUSTIC CHANNEL 


Consider a genera] range-dependent environment such that small regions in the 


vicinity of the source and in the vicinity of the vertical array are assumed to be range 
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independent, but the larger region between is allowed to have arbitrary range dependence 
in bathymetry and sound speed. In the ideal case, the vertical array spans the entire water 
column with elements having uniform spacing d and the modal eigenfunctions have 
negligible amplitude in the bottom. Ignoring absorption losses, the TRA field in the 


vicinity of the source is given by (Kuperman et al., 1998) 


Bai Zi teins Ze) iK,(R)("— R) 


: 2.20 
8p (z,)K,(O)dJ/Rr’ cae 


Prra(’, a @) = S*(wye" S 


where K,,(0) and K,,(R) are the modal eigenvalues evaluated at the vertical array and at 


the source, respectively. ‘V,,(R, z) represent the modal eigenfunctions at the array and d 


is the vertical array inter-element spacing. 

In a situation that is not too far from this ideal case, Eq.(2.20) can be used to illustrate 
the properties that an actual array may possess. In this ideal case, the first property that 
one can notice from Eq.(2.20) is the independence of the focus field with respect to the 
range-dependent environment between the focus and the array. The focus depends only 
on the local properties of the water column and the sea floor, but it 1s not affected by 
bathymetry or range-dependent water column properties in the region between the array 


and the focus. Therefore Eq.(2.20) shows that, neglecting absorption losses and the 


cylindrical spreading factor , the focus pattern is independent of the distance 


—f_ 
VRr’ 
between the source and the array. Thus all the general features of time-reversal acoustics 


discussed in this thesis should also exist in more general, range-dependent environments. 


Losses due to absortion and scattering are expected to cause attenuation of higher- 
order modes. As a consequence, the focus becomes more blurry than that in the previous 
ideal lossless situation. Furthermore, as the distance between the array and the focus 
increases, the focus becomes blurrier due to the strong range and mode number 
dependence of attenuation. However, such influences will not be specifically treated in 


this analysis. 





HI. NUMERICAL MODELING OF THE TIME-REVERSAL ARRAY 


This chapter presents an overview of the parabolic equation (PE) model used for 
acoustic propagation modeling, as well as the method used for its implementation - the 
split-step Fourier (SSF) method (Hardin and Tappert, 1973). The Monterey-Miami 
Parabolic Equation (MMPE) model (Smith, 1996) is the acoustic propagation model used 
throughout this thesis. This propagation model is an upgraded version of the University of 
Miami Parabolic Equation (UMPE) model (Smith and Tappert, 1994). It is based on a 
parabolic equation approximation to the wave equation which was first introduced to the 
underwater acoustics community by Tappert (1974). This approximation is valid for lower 
frequencies (down to 2Hz) than the geometrical acoustic (small wavelength) 
approximation as it retains all the diffraction effects associated with the ocean sound 
channel (Tappert, 1977). 

Other propagation methods such as the separation of variables are based on the 
approximation that the ocean is exactly honzontally stratified. In contrast, the parabolic 
equation model retains full coupling between local waveguide modes, making it valid for 
more realistic, non-stratified oceans. One of the main advantages of the PE/SSF 
implementation is the speed with which the acoustic field can be computed. These are 
some of the reasons why the use of parabolic approximation methods has become 
commonplace in the prediction of underwater acoustic propagation. 

A. PARABOLIC EQUATION THEORY 
The majority of the ocean environment is well suited for a description in cylindrical 


coordinates, particularly in the case of shallow water propagation where the range is much 
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larger than the ocean depth. For this reason, a cylindrical coordinate system is favored 
over others. 


A time-harmonic acoustic field can be described in cylindrical coordinates by 


P(r, z,9,t) = pdr, ee (3.1) 
where pr, z, 9) is the frequency-dependent pressure amplitude as a function of radial 


range r, depth z (positive downward), and azimuthal bearing 9. Substituting P(r, z, 0, t) 
into the inhomogeneous wave equation with a point source leads to the Helmholtz 


equation in cylindrical coordinates (Jensen et al., 1994) 


2 2 
19 PA 4 LOPE SP 


es | ape ee eae 2n2 ee ee 
ror oF Oo? Oz? + kon*(r, 2; O) Pp = 4nP,O(r rs), (3.2) 


Co 


c(r, z, 9) 


where n(r, z, 0)= is the acoustic index of refraction, c, is the reference sound 


7 wm. 
speed, c(r, z,@) 1s the acoustic sound speed, and k,= — is the reference wavenumber. 
Co 


The right-hand term represents a point source at location r, = (7.,2.) = (0, Zz), witht 


pressure amplitude P,, defined at a reference distance lr — r| = Im. Note that c(r, z, 0) 


includes all the features of the environment, although density variations were neglected in 
this derivation. These density variations could be incorporated into a new index of 
refraction without any loss of generality. 

Assuming the ocean acts aS a waveguide, acoustic energy primarily propagates 
outward from the source in the horizontal direction. Therefore, the pressure field can be 


approximated by 


pAr, 2,0) = Wr, z O)H§)(k,r), GB) 
Where Wr, z,@) 1s a slowly varying amplitude function that modulates the outgoing 
zero-th order Hankel function of the first kind, H‘/)(k,r) . Using the far-field asymptotic 


expansion of the Hankel function (Gradshteyn and Ryzhik, 1994), the pressure field 


Fq.(3.3), may be defined by 


R iKor 
PA 2% 9) = p.,. ey Zoe”. (3.4) 


This equation provides the relationship between the acoustic amplitude p,r, z, 0) and the 


so-called “PE field function” Wr, z, 0). The PE field function is normalized such that at 


r=R,, |w4 = 1 and [pj = P,. Neglecting the cylindrical spreading = factor and 


ko 


the range-dependent factor e! ” the PE field function yr, z,@) can be considered as 


being the ocean transfer function at location (r,z,@) due to a source at location 
cs Ly 0) : 
Substituting Eq.(3.4) into the Helmholtz equation and separating the equations for the 


incoming and outgoing fields (e.g., Smith and Tappert, 1994), one can show that the 


defining equation for the PE field function is 


CAVA 2.6) = (1-Qop)UAr, 2 0), 3.5) 


where Q,,,, is the so-called “square-root operator” defined by 


H 


i do 2 2 
QO, = et -1) +1 (3:0) 


Tappert (1977) introduced several approximations to Q,, including the standard parabolic 
equation (SPE). One of the key features of the SPE approximation is the separation of the 


operator Q,,, into two independent operators, 7,, and U 


op a separating the depth 


op 
derivative and environmental term, respectively. The equation defining the evolution of 


the field function y,7, z, >) then takes the form 


2 


EAWAT 0) = (Top 5) + Von) WA 0. 3.7) 


Thomson and Chapman (1983) introduced a higher order approximation based on an 
operator splitting by Feit and Fleck (1978), valid for higher angles of propagation and 


commonly known as the “wide-angle” (WAPE) approximation. In this approximation, 


= i ge i 0 
T,) = ~janea|(2 + aac) + 1| (3.8) 
and 


Ope eee (3.9) 
Jensen et al. (1994) suggests that this approximation extends the accuracy of the solutions 
to a half-bandwidth of 40° but benchmark testing using the WAPE approximation 
demonstrated the capability to accurately ane fields for select environments with 
half-beamwidths greater than 70° (Chin-Bing et al., 1993). The WAPE approximation 


also features less sensitivity to the choice of c, (reference sound speed), and the effects of 
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phase errors are greatly reduced in typical deep ocean conditions relative to the SPE 
approximation (Jensen et al., 1994 and Chin-Bing et al., 1993). For these reasons the 
WAPE approximation is the one used in the MMPE propagation model. 
B. PARABOLIC EQUATION NUMERICAL SOLUTION 

In order to numerically solve the parabolic equation, the MMPE uses a split-step 
Fourier (SSF) algorithm (Hardin and Tappert, 1973). This is a marching algorithm that 


integrates the solution in range by applying the operator U,, in the z-domain and the 


operator T,, in the k,-domain. The latter operator is defined in the k,-domain (vertical 
wavenumber domain) as 
. fe 
Top(k,) = 1- t - : (3.10) 
Note that both operators are just scalar multipliers and may be applied independently. 
The approximate solution for y,(r, z) is iterated in range from rto r + Ar according 


to (Smith and Tappert, 1994) 


—tkoArTop(r, k.) 


wAr+Ar,z) = ge e "(ur 2)))] (3.11) 


where F and F"! represent forward and inverse Fourier transformations, respectively. Note 
that Eq.(3.11) provides time-harmonic solutions to the PE field when substituted into Eqs. 
(3.4) and (3.1). Broadband results are obtained by solving Eq.(3.11) for an evenly-spaced 
discrete set of frequencies in the bandwidth of interest. For each time-harmonic solution, 


the output of the model is in the form of the complex PE field functions for each frequency 


and spatial grid point, y,7;, z;) , referenced to a unit magnitude at r = 1m, and computed 


at the spatial grid points (7;,2;) . 
C. MONTEREY-MIAMI PARABOLIC EQUATION MODELING 

The MMPE acoustic propagation model allows user-specification of the most 
important environmental parameters as inputs. In the model, the environment is divided 
into three distinct layers: the water column, the sediment layer and a deep bottom layer. In 
the water column the user can specify the range-dependent sound speed profile. For the 
sediment layer, the water/sediment bathymetry and several range-dependent acoustical 
parameters can be specified: the sound speed, the sound speed gradient, the density, 
compressional attenuation, shear speed and shear attenuation. All these parameters are 
depth-independent with the exception of the sound speed that has the specified gradient. 
For the deep bottom layer, the user must specify the same parameters as in the sediment 
layer. For a single bottom composition, the deep layer can be defined at a depth below the 
maximum computational depth thereby removing its influence from the calculation. 

The MMPE model allows the use of two types of sources. A wide-angle source which 
approximates a point source and a vertical line array where steering is allowed. Source 
parameters that need to be specified are the depth, the center frequency, the bandwidth, 
and the number of discrete frequencies (should be a power of two for efficient FFT 
computation). The principal input data file specifies the name of the environmental data 
files, the name of the output binary file and other parameters such as the computational 
grid size, the computational range and depth, the output range of values and its grid size, 
and the reference sound speed. The output binary file consists of the PE field function 


evaluated on the output grid at the discrete set of frequencies. It includes a header with the 
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necessary environmental information for the post-processing of the PE field function. For 
more information about the contents of the environmental data files and for the source 
code in Fortran, the reader is referred to the Ocean Acoustics Library supported by the 
U.S. Office of Naval Research through the web page http://oalib.njit.edu/pe.html. From 
this web page the reader may also download the UMPE model documentation (Smith and 
Tappert, 1994). 

D. TIME ARRIVAL STRUCTURE 


It is assumed that the source transmits a single pulse s(t) centered at frequency f, with 


a spectral shape S(f) over bandwidth BW. Note that the bandwidth BW corresponds to the 


bandwidth over which the acoustic field solutions need to be evaluated. The output of the 


MMPE model is the PE field function wr, z) evaluated at N discrete frequencies in the 





bandwidth of interest, from f. — af to f+ af ) , Where f. is the carrier frequency 


and the frequency spacing Af is given by Af = raat Due to the Founer transform 


convention assumed, the time arrival structure at some distance r from the source requires 


the Fourier transformation of the pressure field pr, z) . From Eq.(3.4), the pressure field 


~ is related to the PE field function by 


R ‘kor 
PAY, Zz) = p, fey (r, Ze *e ; (3.12) 
i 


Taking into account the transmitted pulse spectrum, the corresponding pressure field 


becomes 


R ikyr 
PAr, Z) = p, Be (7 Ziseune (sels) 
fi 
The arrival time structure of the propagating field may then be computed from 
_ ve " ikar —i2nft 
Dire. eee = WAr, z)Sfje  e ae (3.14) 
ikr me z 
Noting that we may write e “ =e  °, the result of this phase factor in Eq.(3.14) 
produces arrival time results in terms of “reduced time”’ 
T=t-—, (3.15) 
Co 
where f is the absolute travel time of the pulse. Then 
R, ~i2nfT 
Dee ae [ef War, 2Sfoe ~~ df. (3.16) 
r 


Finally, to decrease the computational burden associated with large transform sizes in 
evaluating the full Founer transform, the MMPE model heterodynes the signal, 


downshifting the center frequency to d.c., and integrates only over the bandwidth such that 


BW 

Z ae P 

i » | Vpag (ts D9 + fer af, 3.17) 
BW 
oF 


where f = f—f, and D(r, z, t) is the envelope of the pressure field. This will introduce a 


Tf 


phase factor of ee" in the time domain which can easily be added later. 
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At this point it is worth mentioning that NV, the number of discrete frequencies and also 
the transform size, should be a power of two in order to use computationally-efficient FFT 


algorithms. Since the signal bandwidth is BW and the transform size is N, then 


jl 
R= ee 
giving a time window of duration 
N-I 
Tyindow = (N—1)- AT = (3.19) 


This time duration value turns out to be very important as a pulse propagates in the ocean 
and “stretches” in time due to multipath propagation. Therefore according to the 
bandwidth, the number of frequencies N should be chosen judiciously in order to avoid the 


undesired “wrap-around” effect due to pulse elongation greater than the time window. In 


practice, the minimum value for N is N = 2”, where L = [log,(BW)]+J and the 


delimiters invoke rounding up to the next integer. 

E. TIME DOMAIN ANIMATION OF SPATIAL ACOUSTIC PROPAGATION 
For a better understanding of the mechanics of time reversal acoustics, it is very useful 

to look at time-domain snapshots of the acoustic field. In order to perform this analysis, 


r 


the “reduced time” at each range r; is converted to absolute time t;, = T,+—, 
Co 


converting p(r;, Z;, T,) to p(7;, Z;, t;,) . The acoustic field can now be considered as being 
evaluated on a time-space grid defined by (7,, Z;, t;,) where t;, represents absolute time. 


Time-domain snapshots are then obtained by defining a specific “time-frame” such that 


t;,=1, and then evaluating p(7;, Z;, t,) at the ranges and depths of interest. Note that since 


il 


the width of the window T,,,, 75, 18 finite, the total range extent R covered at any specific 
Uae Le ey cy vey pares 

To obtain a “time-frame” in the range of interest, some constraints are required on the 
bandwidth of the signal BW and the range step of the output Ar. In general, the time 


difference between two grid points in range is 


Lae = Ui ia orn (3.20) 
The condition for a constant time-frame over different range points is 


laiisayk = lick+m)- (3:21 


Therefore, Eq.(3.20) becomes 


For a given time resolution Ar, one can then have 


Ge - 
Ar = mc,At = Maa Witlii=102.5... 3.233 


where c, is the reference sound speed and m is an integer different from zero. From 


Eq.(3.23), for a given BW, the maximum range resolution allowed in a time-frame is 
C 
= ie For graphical reasons, a Ar that provides about 100 points in the range of 


interest gives good results. 
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IV. ENVIRONMENTAL INFLUENCES IN TIME-REVERSAL ACOUSTICS 


This chapter presents the results of PC array numerical simulations using the MMPE 
acoustic propagation model. Section 4.A gives a brief description of the UWA channel 
used throughout this chapter and Section 4.B gives the description of the transmitted 
signal envelope. Section 4.C gives a description of PC array focusing in the environment 
previously described and how the TRA technique is similar to matched field processing. 
Section 4.D presents a study of the effects that different carrier frequencies have on the 
focus properties. Section 4.E presents the frequency characterization of the channel 
transfer function in order to evaluate the possibility of bandwidth degradation in the 
channel. In Sections 4.F and 4.G several PC array configurations are considered. Section 
4.F shows the effects of changing the element spacing and Section 4.G shows the effects 
of changing the PC array length (or array aperture). Section 4.H presents an 
implementation of focus range shifting. Depending on the source depth, the focus will 
change its location in range when the PC array changes the carrier frequency of the 
transmitted signals. 

During the TRA examination process that lead to this thesis several other 
environments were considered. In all these environments it was observed that the PCA has 
similar focusing properties as the ones described for the UWA channel considered in this 
chapter. For conciseness, only the environment described in Section 4.A will be examined 
here because it is the most realistic of the cases considered. The other UWA channels that 
were considered include channels with different sound speed structures and different 


range-independent as well as range-dependent bathymetry. 
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A. ENVIRONMENT CHARACTERIZATION 

The UWA environment used throughout this chapter to study the properties of PC 
arrays 1s presented in Figure 4.1. This UWA channel is range-independent because the 
sound speed profile, bottom bathymetry and bottom characteristics do not change along 


the range. As depicted in Figure 4.1, the source is located at a range of Om and a depth of 
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Figure 4.1. UWA Channel Profile. 


150m, and the vertical PC array is located 6km from the source. The PC array in Figure 
4.1 represents an illustration of its location with respect to the source. Throughout this 
chapter several array configurations will be used. In some cases the PC array may have 
different element spacing, in other cases it has different lengths and in some cases it can 


span the entire water column. 
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Figure 4.1 also shows the sound speed profile and the bottom characteristics. This 
environment has a upward refracting sound speed profile from a depth of 0m to 50m. 
From 50m to the bottom at 200m, this UWA channel exhibits a sound speed profile similar 
to a typical shallow ocean sound speed structure with the sound axis (minimum value of 
the sound speed) located at 150m in depth. The bottom has compressional attenuation of 
0.1dB/m/kHz, a sound speed gradient of 1s! and a density of 1.5¢/cm?. 

B. TRANSMITTED SIGNAL ENVELOPE 

In most of the TRA numerical simulations presented in this chapter, the transmitted 
pulse has a carrier frequency of 800Hz. The exceptions are in Section 4.D where carrier 
frequencies of 400Hz and 1200Hz are used. In each case, the transmitted pulse has a 
bandwidth of 100Hz with a -3dB effective bandwidth of 36.6Hz. The MMPE 
computational bandwidth over which the solutions of the PE field functions are computed 
includes the full 1OOHz. The envelope characterization of the transmitted pulse is given in 
Figure 4.2. It should be mentioned that the spectrum of the transmitted pulse envelope is 
given by the coefficients of a Hanning window (Proakis and Manolakis, 1996). Therefore 
the lower panel in Figure 4.2 represents a Hanning window (in dB units) that spans 100Hz 
in frequency. The time domain representation of the corresponding signal at a carrier 
frequency of 800Hz can be found in the third panel of Figure 1.2 on page 6. For this 
signal, the -3dB pulse width is 14.1ms. 


C. TIME REVERSAL ACOUSTICS IN A  RANGE-INDEPENDENT 
ENVIRONMENT 


For the environment described in Section 4.A, Figure 4.3 represents the magnitude of 
the ocean response for a time-harmonic acoustic source at a frequency of 800Hz in a 


vertical plane defined by the source and the array elements. The ocean frequency response 
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Figure 4.3. Ocean Response at 800Hz. 
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at 800Hz is represented in terms of transmission loss (dB re 1m). Considering the ocean as 
a spatial filter, Figure 4.3 represents the “forward propagation transfer function” at the 
frequency of 800Hz. The term “forward propagation” is used to denote the propagation 
from the acoustic source. If the source located at a depth of 150m transmits a pulse at the 
carrier frequency of 800Hz with a bandwidth of 1O0QHz (-3dB effective bandwidth of 
36.6Hz), the time arrival structure at the vertical array may be evaluated by Fourier 
synthesis of single frequency (CW) solutions over the bandwidth of interest. The time and 
frequency domain characterization of the transmitted pulse was previously described in 
Figure 4.2. 

Figure 4.4 represents the time arrival structure of the received signal envelope at the 
PCA for the 800Hz signal with 100Hz bandwidth. For a given array element, the envelope 


of the received signal is given by the corresponding horizontal line at the array element 
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Figure 4.4. Time Arrival Structure at the PCA. 


on 


depth. The effects of multipath propagation are clearly observed in Figure 4.4 where any 
array element detects several multipath arrivals. From Figure 4.4 one can see that the 
multipath structure 1s similar between closely spaced array elements. The diversity 
between array elements can be exhibited by computing the envelope cross-correlation of 
the time arrival structure represented in Figure 4.4. The envelope cross-correlation shown 
in Figure 4.5 was computed with respect to the envelope received by the array element at a 


depth of 150m. With respect to this array element one can observe that the other array 
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Figure 4.5. Envelope Cross-Correlation across the Array Elements. 
elements become more decorrelated as the depth difference between array elements 
increases. In general as the distance between array elements increases the multipath 
structure becomes more decorrelated. This spatial decorrelation between array elements is 
due to the presence of multipaths and is used in communication systems that explore the 


spatial diversity in the water column. 
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At the PCA, the received signals are time-reversed and transmitted back to the source. 
Figure 4.6 represents the ocean frequency response to the time-reversed acoustic field 


transmitted by the PCA evaluated at 800Hz. At this point it should be mentioned that the 





Transmission Loss (dB re 1m) at Freq = 800Hz. 
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Figure 4.6. Backward Propagation Ocean Response at 800Hz. 
term “backward propagation” will be used throughout this thesis and by that is meant 
propagation of the time-reversed signals from the PCA through the reciprocal 
environment. From this perspective, Figure 4.6 represents the back propagation ocean 
frequency response at 800Hz. 
In Figure 4.6, the propagation is from the left to the right, i.e. the PCA 1s located at 
range Okm and spans the entire water column. In this geometry the original source is 
located at a range of 6km and at a depth of 150m. As expected, the phase conjugation of 


the 800Hz frequency component and its consequent propagation through the channel 


produces a strong spatial focus at the source location. 
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From the viewpoint of linear systems theory, the ocean itself behaves similarly to a 
spatial matched filter. Due to spatial reciprocity, the only location in space where this 
spatial filter is matched to the forward propagation function is at the source location. 
Therefore the output of this spatial filter is maximized at the source location creating a 
strong focusing effect because the source location is the spatial location where the 
backward propagation transfer function matches the forward propagation transfer 
function. In other words, at the source location all multipath contributions will arrive in 
phase and add constructively creating the focusing effect. Using these same arguments, all 
the other frequency components of the transmitted pulse will focus at the source location. 


Figure 4.7 represents the time arrival structure of the back propagated pulse (f.=SO0Hz, 


BW=100Hz) at the source range and the temporal, vertical and horizontal PCA focusing 
properties. The data displayed corresponds to the recorded signals of a vertical array of 
elements that spans the entire water column at the source range. 

In Figure 4.7 one can observe the temporal focusing at the source location where the 
multipath structure is greatly suppressed. This temporal focusing is also a consequence of 
the matched filter behavior of the ocean where all frequency components will arrive at the 
same time at the same location. In this numerical simulation the PCA spans the entire 
water column with an element spacing of d = 0.244m and A/d = 7.67, approaching the 
ideal case of a continuous line of sources. As a consequence, the temporal sidelobes 
represent the limiting situation that can be achieved in practice. As will be shown in this 
chapter, decreasing the number of array elements increases the temporal sidelobes. 

Also note in Figure 4.7 that this focusing process has vertical (along the depth) 


sidelobes. As mentioned in Chapter IH, the vertical focusing is a consequence of the 
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closure property of the propagating modes. Thus, these sidelobes appear in part because 
there are higher order modes, called evanescent modes (Jensen et al., 1994), that decay 


rapidly in range and don’t propagate in the channel. In addition, some propagating modes 
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Figure 4.7. Time Arrival Structure Focus Range. Temporal, Vertical and Horizontal PCA 
Focusing Properties. 


are attenuated differently by the environment. Thus the focusing extent is dependent on 
the number of modes allowed in the channel. 

Also displayed in Figure 4.7 is the horizontal focusing at the source depth. This plot is 
very representative because it gives the maximum amplitude of the received signal 
envelope for a given location in range. The main lobe is located at 6km corresponding to 
the source range and the -3dB main lobe width is about 60m, much larger than the 


corresponding value along the depth which is 3.3m. 
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D. FOCUSING AT DIFFERENT FREQUENCIES 


In this case study, three signals with carrier frequencies of 400Hz, 800Hz and 1200Hz 


were used. All these signals have the same bandwidth of 100Hz that corresponds to a -3dB 


effective bandwidth of 36.6Hz, the same envelope spectral shape as described by Figure 


4.2 and the PC array spans the entire water column. The temporal, vertical, and horizontal 


focusing properties for each signal are displayed in Figure 4.8. 
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Figure 4.8. Focusing at Different Frequencies. 


As shown in Figure 4.8 the temporal duration of the main arrival is identical for any of 


these test cases. This is a significant result because it means that using the same bandwidth 


at different carrier frequencies achieves the same temporal pulse resolution and 


consequently the same symbol rates. It appears that the channel frequency response does 


42 


not change significantly over the bandwidth of each signal and therefore the same 
temporal resolution exists at these different carrier frequencies. The temporal sidelobes 
are below 40dB and decrease with increasing frequency. This is presumably due to the 
increase in the number of propagating modes at higher frequencies. In addition, the focus 
dimensions (or “footprint’) are found to decrease with increasing frequency. This is due to 
a combination of the increase in the number of propagating modes and the decrease in 
acoustic wavelength. 
E. CHANNEL CHARACTERIZATION IN THE FREQUENCY DOMAIN 

For the channel described in Figure 4.1 where the source is 150m deep and the PC 
array spans the entire water column, Figure 4.9 represents the forward propagation 


channel frequency response evaluated at the PC array element 150m deep. The upper 









aor Channel frequency response at PCA (depth=150m). 
7 [7] | a a = aoe ea TA [== a aan aa 
1.6 
14 
12 
o 
a f | 
2os 
= i | : 
06 | | | 
04 a ‘il HN ll 
0.2 
0 ee Si ed = 5 —— 
200 400 600 800 ~<a 000 1200 





Frequency (Hz) 


Channel pene cy FesPOnse i at PCA (cen ae (in dB). 
-130;-—— ia =e a 








=e =400Hz _ 
f =800Hz 
_ f =1200Hz 


eE 4 oe | 
-20 -10 0 10 20 30 50 
Frequency in respect to each carrier frequency (Hz) 








Figure 4.9. Forward Propagation Channel Frequency Response at the PC 
Array Element Located 150m Deep. 
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panel represents the frequency response from 200Hz to 1400Hz and the lower panel 
represents the channel frequency response in the range of frequencies used in the previous 
section. Regarding the lower panel one can observe that the channel frequency response 
has smooth peaks and deep fades almost evenly spaced. In some cases the deep fades can 
be 50dB below the peak level. This behavior in the forward propagation frequency 
response 1s a consequence of the multipath interference structure at this array element. For 
the other PC array elements the channel frequency response exhibits similar behavior with 
the deep fades located at different frequencies that the ones shown in this example. 

Figure 4.10 represents the channel frequency response at the focus location. This 
channel frequency response will be designated by “two-way” frequency response because 
it includes forward propagation followed by backward propagation to the focus location. 


Comparing Figure 4.9 with Figure 4.10, the two-way propagation frequency response 
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Figure 4.10. TwoWay Channel Frequency Response at Focus Location. 
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evaluated at the focus location exhibits a smooth behavior with frequency where the deep 
fades are absent indicating that all frequency components at the focus location are in 
phase. Note however that the lower frequency components at the focus have lower 
amplitude than higher frequency components. This is presumably due to stronger 
interaction of the low frequency modes with the lossy bottom. 

In Figure 4.9, only the frequency response for the 150m deep array element is shown. 
Since the PC array spans the entire water column and the interference structure will 
change along the depth, the other PC array elements will exhibit nulls in the frequency 
response at frequencies other than the ones shown for this PCA element. When back 
propagation is performed, the contributions from all the PC array elements will add 
coherently at the focus location, restoring all frequency components in the received field. 
Thus, at the focus location the two-way frequency response is much smoother than for the 
forward propagation case. 

FL EFFECT OF INTER-ELEMENT SPACING 

In this section we evaluate the influence on the focus due to PC arrays with different 
element spacing. In this test case the transmitted signal has a carrier frequency of 800Hz 
and a bandwidth of 100Hz. The geometry of the different PC arrays used in this study is 
given in Figure 4.11 where each array is identified by a number from | to 6. Associated 
with each PCA are 3 numbers giving the number of elements, the distance between 
consecutive elements in meters and the distance between consecutive elements in terms of 
carrier wavelengths. Note that PCA number | has the smallest element spacing and that 


the distance between array elements doubles between consecutive arrays. 
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Figure 4.11. Focusing Properties Using PC Arrays with Different Element Spacing. 
Starting from PCA number | to number 6, as the element spacing increases, the array 
cannot provide enough spatial sampling of higher-order propagating modes. Thus, the 
focusing degrades as a result of the lack of resolution of higher order modes in the field. 
From a linear systems point of view the changing in focus is the result of removing array 
elements that are more correlated, i.e. as the array spacing increases it increases the 
decorrelation between array elements. 
From Figure 4.11 one can conclude that the spatial dimensions of the focus do not 
change significantly for these PCA configurations. Comparing these results to Figure 4.7 


where the PCA spans the entire water column one should expect an increase in temporal 
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sidelobes because the PC arrays in Figure 4.11 only span half of the water column, 
collecting less information about the channel conditions. Remarkably, in the worst case 
situation (PCA number 6 with only 9 elements), the temporal sidelobes are 25dB below 
the main lobe. These results seem to indicate that even with a small number of array 
elements, a good focusing quality can be achieved provided that the PCA has enough 
spatial sampling of the propagating modes in the water column. 
G EFFECT OF APERTURE SIZE 

In this study the element spacing is kept constant in order to evaluate the effect that 
different PCA aperture sizes have on the focus quality. The geometry and spatial location 
in the water column of the different PC arrays is depicted in Figure 4.12. In this figure, 
each array is designated by a number from | to 5. The numbers associated to each array 
give the number of array elements and the array length in meters and in carrier 
wavelengths. As the PCA aperture size increases, it can provide a better sampling of the 
low-order propagating modes. Starting with array number | and going to array number 5, 
the changes in focus are due to better spatial sampling of the low-order propagating 
modes. Thus, one should expect better focusing properties as the aperture size increases 
because the PCA has more spatial diversity and provides a better sampling of the channel 
conditions. 

Because these PCA arrays cannot properly sample low-order propagating modes, the 
sidelobes increase in the depth and range coordinates. In Figure 4.12 the first three arrays 


5X. and 10X.., respectively, have temporal sidelobes about 12dB 


with apertures of 2.5A.., 5 
below the main lobe. When the aperture size increases to 20A,. and 40A.., the temporal 


sidelobes decrease to 27dB and 35dB, respectively, below the main lobe. This distinctive 
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Figure 4.12. Focusing Properties Using PC Arrays with Different Aperture Size. 


behavior of array numbers 4 and 5 compared to the other three is because they provide a 
better sampling of the lowest-order propagating modes and therefore have better temporal 
properties. 

Presumably, the placement of these arrays in depth will also affect the focus quality. 
The best quality focus would be expected when the array is placed at a depth which 
optimally samples the dominant propagating modes. The conclusion that can be drawn is 
that in a given environment one should evaluate the most important propagating modes 


and use an array aperture that provides enough sampling of these modes. 
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H. EFFECT OF FREQUENCY BIN SHIFTING ON FOCUSING AT OTHER 
LOCATIONS IN RANGE 


The basic idea of the range shifting approach is that the PCA can focus at different 
ranges by increasing or decreasing the carrier frequency of the transmitted signal. This 
frequency shifting of the carrer frequency can be implemented in near real time by an 
FFT bin shift at the PCA array before transmission. 

This procedure 1s based on the acoustic-field invariant property in the coordinates of 
frequency and range (Brekhovskikh and Lysanov, 1991). [t has been shown (Chuprov and 
Mal’tsev, 1981, and Chuprov, 1982) that the interference structure in oceanic waveguides 
is characterized by the existence of lines of maximum intensity /(@, r) having a slope B, 


and this slope 1s an invariant for a particular group of modes (Song et al., 1998), 





daa eee CLAD 
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where v and uw are the phase and group velocities characterizing a given mode group, 
respectively. The invariant B is determined by the properties of the channel and the 
derivative in Eq.(4.1) averaged over this group. 

This invariant property applies not only to the sound intensity but also to the phase 
velocity and envelope (group) delay (Grachev, 1993) thus extending to broadband pulse 


propagation. From Eq.(4.1), if the PCA transmits the complex conjugated field at a carrier 


frequency of @. + Ao), the corresponding range shift is given by 


af 
(om — AC (4.2) 
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The results in Figure 4.13 show that range shifting was attained when the PCA 
receives a signal at a carrier frequency of 800Hz with 100Hz bandwidth and shifts the 


carrier frequency up or down by 2.5Hz and 5Hz. These results show that when the carrier 
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Figure 4.13. Focus Range Shifting by Changing the Carrier Frequency at PCA. 


frequency is shifted by +2.5Hz, the focus range shifts 30m with respect to the focus 
location at the nominal carrier frequency. When the carrier frequency is shifted by +5 Hz, 
the focus range changes by 60m. These results give an approximate value for the invariant 
8B of 0.119. This range shifting property was also observed when the source was located 
50m and 100m deep and by shifting the carrier frequency by +2.5 Hz and +5 Hz, giving a 


value for B of 0.228 and 0.258, respectively. The results in Figures 4.14 and 4.15 


SO 


correspond to a source depth of 50m and 100m, respectively. At this point it seems clear 


that this invariant B has depth dependency. 


It should be mentioned that the range shifting property can be observed as long as Aw 


is small when compared to the transmitted signal bandwidth. In other situations where 


+20 Hz and +40 Hz shifts were used, some dissimilar results were found. For these carrier 


shifts, no range shifting was observed when the source was Jocated 150m deep. When the 


source was located 50m and 100m deep, it was observed that the focus split into two foci. 


These two foci are symmetrical] in range with respect to the nominal focus range (no range 


shifting) and are also symmetrical in depth with respect to the nominal! focus depth. 
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Figure 4.14. Focus Range Shifting for a Source 50m Deep. 
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Figure 4.15. Focus Range Shifting for a Source 100m Deep. 
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V. APPLICATION OF TIME REVERSAL ACOUSTICS TO 
SHALLOW WATER COMMUNICATION SYSTEMS 


This chapter illustrates the application of TRA to UWA communication systems. It 
illustrates the use of PC arrays both in a point-to-point acoustic link and in a network 
situation. The communication system described in this chapter uses a reduced-complexity 
non-coherent detector at the expense of increased complexity in the transmitter (PCA). 
The examples illustrated in this chapter correspond to a nominal signaling rate of 30.5 
symbols/second where 4 bits of information are encoded in each symbol, as described in 
Section 5.C. The channel used in this chapter for numerical simulations is the one 
described in Section 4.A. 


A. USING PHASE-CONJUGATED ARRAYS IN UNDERWATER ACOUSTIC 
COMMUNICATIONS 


In this section, an example of the use of PC arrays in UWA communications Is given. 
Figure 5.1 describes the simplest case of point-to-point communications between 2 
stationary communication nodes. Each communications node of Figure 5.1 has a PC array 
that spans the entire water column and the source/receiver element of each node is located 
at a depth of 150m. 

As seen in the previous chapter, when the time-reversed acoustic field propagates in 
the channel, it will focus back at the original source location due to spatial reciprocity. 
From the viewpoint of systems theory, the channel can be considered a spatial filter. In the 
forward propagation case (acoustic propagation from a source to a PCA), the presence of 
multipath and other channel distortions causes a different channel impulse response at 


different locations in the channel. After back propagation from each element of the PCA, 


2%, 


the channel behaves like a matched filter where the original source location is the only 
location in space where the forward propagation transfer function matches the backward 
propagation transfer function for each PC array element. Therefore during 
backpropagation the output of this spatial filter will be maximized at the source location. 
From these observations one can conclude that the source position is the optimal location 
in space to place a receiver. 

For each node in Figure 5.1, the PCA focusing properties require the acoustic source 
to be co-located with the designated receiver element for TRA UWA communications. 
This receiver element will be responsible for reception of the TRA UWA communication 
signals. At this point is possible to conceive that each node will be using a single-element 


transducer that can be operated both as an acoustic source and as an acoustic receiver. 
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Figure 5.1. Channel Geometry in a Point to Point Communication Link. 
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The first step towards the focusing of PCA 1 requires a probe (pilot) pulse 
transmission from the node 2 acoustic source. As this probe pulse propagates, it will 
incorporate channel distortions including multipath. In the next step, the received acoustic 
signals at each PCA | element are digitized and stored in digital memories. For PCA 1 
transmission (or backpropagation), the recorded signals will be time-reversed and 
transmitted by each PCA | element. The time-reversal transformation at PCA 1 is a simple 
task to perform, requiring the digital memories to be read in a last-in-first-out (LIFO) 
manner. At the PCA transmission, each element is excited by the output of the 
corresponding memory where the received signals were recorded. 

In order to attain good focusing quality there is the requirement that all PCA elements 
must be synchronized at reception and transmission. By synchronization is meant that it 
will record the received signals during the same time interval and at transmission all the 
elements are fired at same time. 

At this point one can consider that PCA 1 is “focused” at the source/receiver location 
of node 2 because it has recorded the necessary environmental information to generate a 
spatial focus at the node 2 source/receiver location. A similar process can be described for 
the focusing of PCA 2 that requires a probe pulse transmission from the node 1 acoustic 
source. 

The PCA’s focusing procedure is a critical aspect in a TRA communication link. TRA 
communications can only be established after both arrays are focused. In a real situation 
one must take into account the temporal variability of channel conditions. Focusing will 
degrade with time, requiring periodic PCA re-focusing. In a real situation where generally 


the nodes may have some motion, the PC array may need to update the environmental 
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information every time the receiver moves out of focus. In such a siti the range 
shifting property of the PC array discussed in the previous chapter may be a very useful 
feature in a UWA communication system. If the receiver is moving slowly, the PCA is 
able to keep the focus at the receiver location, decreasing the rate at which the UWA 
communication system updates the environmental information (probing). 

As a final remark the electronic circuitry associated with each PCA element includes a 
receiver amplifier, an A/D converter, digital filters, memory, a D/A converter and a 
transmitter amplifier. Until this point no peterenee was made to the presence of digital 
filters at each PCA array element. Its presence will become clear in the next section. 

B. SIGNALING SCHEME 

In this section, a signaling scheme appropriate for PCA non-coherent communications 
is described. _ this signaling scheme each symbol represents a 4-bit word giving a total of 
16 different words in the code. Since this communication system is based on non-coherent 
detection, the value of each bit within a symbol is given by the detection or lack of 
detection of the corresponding frequency component. There are 15 different symbols in 
the code having at least one spectral component. The “all-zeros” word corresponds to the 
case where no carrier frequencies are detected. This signaling scheme corresponds to a 
QFSK modulation scheme. Thus, a simple receiver can be built based on a four element 
filter bank followed by an envelope detector. The detection of a given frequency 
component is achieved when the corresponding envelope detector output is higher than a 
given threshold. 

In this signaling scheme, four frequency components are defined centered at 750Hz, 


800Hz, 850Hz and 900Hz. Each one of these frequency components spans 100Hz in 
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frequency and has a -3dB bandwidth of 36.6Hz. The least significant bit corresponds to 
the frequency component centered at 900Hz and the most significant bit corresponds to 
the component at 750Hz. 

To initialize the PC array, the source at each node must transmit each 100Hz 
bandwidth probe pulse separately. An alternative approach to simplify probe pulse 
transmissions would require that each PCA element has a bank of 4 digital filters centered 
at 750Hz, 800Hz, 850Hz and 900Hz. This way the array can store the frequency 
component contents from the arrival of a single probe pulse covering the entire bandwidth, 
700Hz to 950Hz. 

If the first approach is chosen there is no requirement for a digital filter bank in each 
PCA element but each node must make sure that the probe pulse transmissions will not 
overlap in time at the PCA. 

The choice of which approach to use may be dictated by the temporal variability of the 
channel conditions that will degrade focusing, requiring periodic update of the channel 
conditions. The filter bank approach has the advantage that it requires only one probe 
pulse transmission from each source in Figure 5.1. Thus the communication system 
spends less time in array focusing and higher data throughput can be attained at the 
expense of a more complex PCA. 

C. PHASE-CONJUGATED ARRAY 

In this section it 1s assumed that the PCA is able to learn the environmental 
information in the bandwidth of interest. The PCA is not in a shadow zone and each 
source transmits a probe pulse with bandwidth wide enough to accommodate the 


bandwidth used to transmit all individual frequency components of the signaling scheme. 


Di 


The signaling scheme described in the previous section requires a four element filter 
bank for each PCA element. This four element filter bank frequency response is given in 
Figure 5.2. All the filters have the same spectral shape but different center frequency, each 
filter frequency response spanning over 100Hz giving a -3dB bandwidth of 36.6Hz. The 
frequency response of each filter is specified by the coefficients of a Hanning window. In 
Figure 5.2 the center frequency of each filter corresponds to the center frequency of one of 
the frequency components used in the signaling scheme. 
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Figure 5.2. PCA Element Filter Bank Frequency Response. 

After PCA reception and filtering, each PCA element will have recorded 4 different 
signals with different center frequencies. Each one of these signals will have the transfer 
function of the channel imprinted on it. As seen in the previous chapter, if all PCA 
elements transmit only the signal component centered at 800Hz, a pulse with 800Hz 


carrier frequency and a -3dB bandwidth of 36.6Hz will be received at focus location. 
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D. POINT TO POINT MESSAGE TRANSMISSION 

This section presents a numerical simulation where a message is sent from node | to 
node 2 of Figure 5.1. At each PCA | element, the transmitted signal corresponding to a 
given symbol can be obtained by time-domain superposition (i.e., simultaneous 
transmission) of its different frequency components. Similarly the signal corresponding to 
a sequence of symbols can be obtained by superposition of the signals corresponding to 
the different symbols, properly delayed in the time domain. This time-domain delay 
corresponds to the signaling period and should be greater than the temporal duration of 
each symbol at the receiver location to avoid inter-symbol interference. In this case each 
symbol component will have at the receiver a -3dB temporal width of 14.lms and a 
conservative value of 32.8ms was chosen for the signaling period. This provides a 
signaling rate of 30.5 symbols/second. This choice for the signaling period is somewhat 
arbitrary and based simply on producing results with clearly distinguishable arrivals. With 
accurate processing techniques, smaller signaling periods could be used thereby 
increasing the signaling rate. 

Figure 5.3 presents the spectrogram of the signal corresponding to a message 


composed by the sequence of words from (0001), to (1111), that 1s intended to be 


transmitted from node | to the node 2 receiver. Figure 5.3 represents the ideal signal that 
node | wants to focus at the node 2 receiver. Note that the least significant bit corresponds 
to the frequency component centered at 900Hz and the most significant bit corresponds to 
the component at 750Hz. 

Figure 5.4 presents the results from a numerical simulation where node | spans the 


entire water column and transmits the appropriate signals to focus the message of Figure 
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- 


5.3 at the node 2 receiver. All the signals in Figure 5.4 are collected at a depth of 150m. 
The upper-left panel represents the signal transmitted from the PCA | element located at 


150m depth while the other three correspond to the signal spectrogram at various ranges. 
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Figure 5.3. Spectrogram of the Signal Containing the Desired Message if 
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Figure 5.4. Message Signal Spectrogram Received by a Single Element Receiver 
150m Deep, Evaluated at the PCA (Upper Left), at Two Intermediate Ranges and at 
the Focus (Lower Right). 
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Only the lower nght panel in Figure. 5.4 corresponds to the signal spectrogram at the 
focus location. 

In Figure 5.4 one can clearly observe that as the pulse back-propagates in the channel 
the multipath structure becomes less significant as the range approaches the focusing 
range. It can also be seen that when the pulse approaches the focusing range it increases 
the power spectral density of its spectral components. This increase in the power spectral 
density is a consequence of the vertical, horizontal and temporal PCA focusing properties 
as the sequence of symbols approaches the focus location. 


E. PHASE-CONJUGATED ARRAY VERTICAL DIVERSITY IN 
UNDERWATER COMMUNICATIONS 


Presented in this section is an extension of the TRA technique described in the 


previous sections. The setting for this TRA experiment is depicted in Figure 5.5 and is 
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Figure 5.5. PCA Geometry for Vertical Diversity Examination. 
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similar to the one presented in Figure 5.1 with the difference that node 2 has now 3 source/ 
receiver elements. The extension of the TRA technique presented in this section exploits 
node | PCA vertical diversity to focus the back propagated acoustic field at the different 
node 2 receiver locations. As a result, higher data throughputs can be achieved in the data 
link from node 1 to node 2 by simultaneously transmitting symbols designed to focus at 
different node 2 receivers. 

In this case the PCA 1 focusing process requires probe pulse transmissions from each 
source/receiver element of node 2. The transmission of probe pulses from node 2 has the 
constraint that it must allow enough time between consecutive transmissions so each pulse 
will arrive at PCA 1 at different time slots and the pulses do not interfere. In this example 
PCA 1 will require a data bank with 3 times the storage capacity of the one used in PCA 2 
because it has to store the signals received from 3 probe pulses. The TRA acoustic link 
from node | to node 2 can accommodate 3 times the data throughput of the case presented 
in the previous sections. In this case node 2 has 3 receivers and node 1 is able to send 3 
different symbols (one for each node 2 receiver) during each symbol period. 

At transmission, PCA 1 will assemble the message intended for each receiver by the 
superposition of its frequency components stored in its data bank. The signal transmitted 
by each PCA 1 element will be composed of the superposition of the signals for the 
different node 2 receivers. Figure 5.6 represents the desired message signal spectrogram 
that node 1 wants to focus at each node 2 receiver. Note that these signals do not 
correspond to the actual signals sent by PCA 1 elements. At each PCA 1 element the 


transmitted signal will include multipath and other channel distortions. 
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Figure 5.6. Desired Message Signals for Node 2 Receivers. 


Figure 5.7 presents numerical simulation results of the PCA 1 transmitted signal 
propagation in which Figure 5.6 represents the message that PCA | wants to focus at each 
receiver location. In this simulation, PCA | spans the entire water column and the upper- 
left panel of Figure 5.7 represents the actual signal sent by the element located at a depth 
of 150m. We can observe multipath structure and strong interference between the different 
transmitted messages. The remaining panels in Figure 5.7 represent the signal received by 
each one of the node 2 receivers. Due to the vertical focusing properties of PCA 1, these 
results show very little message intersymbol interference at each receiver and the desired 
message can be easily decoded. 

With respect to this same numerical simulation, Figure 5.8 represents the signal 


spectrogram at different ranges between node | and node 2 collected at a depth of 150m. 
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The results in this figure show that it will be extremely difficult to decode the original 
messages at these locations because the different message signals are overlapping both in 
time and in frequency. These results also show that this TRA vertical diversity technique 
may have applications to UWA secure communication links. Channel distortions and 
signal overlap both in time and frequency of the different messages seems to be a very 
good way to encode digital communication signals with the advantages of a higher data 
throughput. 


F. PHASE-CONJUGATED ARRAY SPATIAL DIVERSITY IN UNDERWATER 
COMMUNICATIONS 


In this section the concept described in the previous section is further extended to 
consider the case of an UWA communication network. The spatial arrangement of the 
acoustic network illustrated in this section can be observed in Figure 5.9. In this case all 
: Channel geometry - spatial diversity 
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Figure 5.9. PCA 1 Geometry for Examination of Vertical and Horizontal 
Spatial Diversity in a UWA Communication Network. 
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the acoustic nodes are located in the same vertical plane. This constraint is only due to the 
fact that MMPE is a 2-dimensional propagation model, and therefore all the sources and 
receivers must be in the same vertical plane. In a real situation one could have these 
acoustic nodes located in different azimuthal directions from node 1. In Figure. 5.9, nodes 
2,3, 4 and 5 are located 2km, 5km, 6km and 7km, respectively, away from node 1. At each 
node the receiver may be located at different depths and node 4 has 2 different receiver 
elements. 

The objective in this section is to illustrate the case where node 1 transmits 
simultaneously in time and frequency different message signals to each one of the 
receivers at the other nodes. In order to accomplish this, node 1 takes advantage of the 
spatial focusing diversity of PCA 1. In other words, PCA 1 is able to simultaneously focus 
different messages at different locations. 

In a network environment, this technique can be difficult to implement. Each PCA has 
to record and store the received probe pulse from each one of the source/receiver elements 
at the other communication nodes. At each focusing stage only one probe pulse is sent by 
one of the source/receiver elements and will be recorded by the PCA’s at the other 
communication nodes. At the end of the probing process PCA 1 will have stored the 
necessary signals to focus at all the receivers. For example, if PCA 1 wants to send a 
message to source/receiver nr. 5, it uses only the stored signals from source/receiver nr.5 
probe signal to compose signals to be transmitted that will focus the desired message at 
the receiver 5 location. In the case where node 1 broadcasts a message for all receivers, 
PCA 1 will use all the signals recorded from each source/receiver probe pulse 


transmission to compose the signals that will be transmitted by each element. 
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Figure 5.10 represents the desired message signal spectrograms that node | wants to 
focus at each receiver of Figure5.9. Note that these signals do not correspond to the actual 
signals sent by the PCA 1 elements. At each PCA 1 element the transmitted signal will 
include multipath and other channel distortions due to spatial variability of channel 
conditions. As seen in Figure. 5.10, four of the desired received signals only have one 


frequency component. The motivation behind this choice of message signals is to observe 
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Figure 5.10. Desired Signals to Focus at the Receivers. 
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intersymbol interference due to temporal overlapping of the different transmitted 
messages. 

Figure 5.11 represents numerical simulation results corresponding to the situation 
where node | transmits simultaneously the appropriate signals to focus each message of 
Figure 5.10 at the corresponding receiver. Again, in this simulation it was assumed that the 
node | PCA spans the entire water column. The results in Figure 5.11 show that little 


interference was attained and each receiver is able to decode the correct message. 
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Figure 5.11. Received Signal Spectrogram at the Receivers. 
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Figure 5.12 represents the received signal spectrogram at different ranges between 
node | and node 2 collected at a depth of 150m. The results in this figure show that it will 
be extremely difficult to decode the original messages at these locations because the 


different message signals are overlapping both in time and frequency. 
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Figure 5.12. Received Signal Spectrogram at a Depth of 150m and Different 
Ranges from Node I. 
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VI. CONCLUSIONS 


In this thesis an attempt was made to characterize the influences of a phase- 
conjugation array design on its focusing properties. Particular attention was given to PCA 
applications in UWA communication systems. This thesis provides a starting point for a 
PC array design, and gives support to the planning of further TRA UWA communication 
experiments both in the laboratory and in an oceanic environment. In the early stages of a 
TRA UWA communication system design, TRA numerical modeling is very important 
because it gives insight to the propagation issues that are involved without the need to go 
out to the ocean and perform very expensive experiments. However, although numerical 
modeling is an excellent analysis tool, it does not provide complete information about 
time-variability of the channel conditions. To take that into account and incorporate it in 
the UWA communication system design, oceanic experiments are planned for the near 
future. 

The numerical analysis of the PCA focusing properties showed that the focus footprint 
decreases in dimensions as the carrier frequency increases. Furthermore, the horizontal 
dimension is larger than the vertical dimension of the footprint. Therefore vertical motion 
of the receiver is much more cmitical than horizontal motion because the receiver can 
easily move out-of-focus. 

Using the same bandwidth and changing carrier frequency in the range 400Hz to 
1200Hz, numerical modeling showed that temporal focusing does not change 
significantly. This is the case illustrated in Section 4.D with an effective pulse bandwidth 
of 36.6Hz. In a more general situation where a larger pulse bandwidth is used, the results 


in Figure 4.10 suggest that for carrier frequencies higher than 1200Hz, some degradation 
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will occur in the pulse temporal resolution at the focus location since the two-way channel 
frequency response will change considerably over the bandwidth. 

Using PC arrays with different element spacing or with different aperture size, the 
PCA footprint did not appear to change its dimensions. Keeping constant the aperture size 
and changing the PCA element spacing, the temporal sidelobes remained approximately at 
the same level. As illustrated in Figure 4.11, temporal sidelobes 25dB below the main lobe 
can be attained with an element spacing of 6.6 wavelengths (in this case the aperture size 
is 100m, half of the channel depth). 

Keeping constant the PCA element spacing and changing the aperture size, the vertical 
and horizontal footprint size did not appear to change significantly. A small improvement 


in horizontal sidelobes was observed in Figure 4.12 as the aperture length was increased. 


Figure 4.12 also showed that with PCA aperture sizes of 2.54, (7 transducers), 5A. (13 


transducers) and 10A,, (25 transducers), the temporal sidelobes remained at the same level 


with respect to the main lobe. A dramatic improvement in temporal sidelobes was 


observed when the aperture size increased to 20A,. (49 transducers) because the PCA was 


able to provide better spatial sampling of the lowest propagation modes. In other words, 
there was more spatial diversity in the PC array. 

What can be concluded from these experiments is that aperture size plays a more 
significant role than the number of array elements in a PC array design. In the examples 
shown in Figures 4.11 and 4.12, the PC array with an aperture size of 100m and only 9 
transducers had lower temporal sidelobes than the PC array with an aperture size of 37.5m _ 
and 49 transducers. The number of PCA transducers is a critical factor in PC array design 


because it will drive up the cost. Each PCA transducer has its own electronic circuitry 
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composed of receiver amplifiers, A/D and D/A converters, storage memory and 
transmitter. In general, underwater devices rely on batteries to operate their electronic 
systems and have severe constraints on the amount of power that can be used. Thus, for 
cost effectiveness one would be better off with a small number of array elements extended 
over a large aperture size. 

Another interesting TRA property that was observed is the horizontal shifting of the 
footprint location when the PCA shifts the carrier frequency of the transmitted signal. At 
this point it seems that this property may be useful in a UWA communication system 
where the stations may be in motion or drifting. Its implementation requires further 
research because each PCA must sense the rate of relative motion in order to keep the 
focus at the receiver location. The relationship between the frequency shifting and range 
of focus shifting is also environmentally dependent. 

For UWA communication systems using noncoherent detection, the use of PC arrays 
seems promising because there is no requirement of carrier phase tracking. The temporal 
focusing properties of the PCA overcomes the requirement of guard times between 
consecutive symbols to ensure that all the reverberation will vanish before each 
subsequent symbol is received avoiding intersymbol interference at the receiver. Thus 
some data rate improvement is made because the insertion of idle periods of time results in 
a reduction of the available data throughput. Array processing is no longer required 
because the PCA footprint is well defined at the receiver location. 

As seen in Chapter V, the data throughput can be further improved by taking 
advantage of the PCA spatial diversity focusing properties. In the cases presented, the 
PCA was able to focus simultaneously different messages at different receiver locations. 


The message signal designated for a given receiver only focused at the intended receiver 
if 


location. At the receiver the corresponding message signal was compressed in time and 
space and had higher magnitude than the other interfering message signals. As a 
consequence, each message signal could only be decoded at the desired receiver. 
Furthermore, this seems to be a good technique for secure UWA communications because 
the different message signals will be overlapping both in time and frequency at other 
locations, making it difficult to recover any of the transmitted messages. 

Due to the PC array spatial diversity focusing properties, PC arrays may have an 
important role in an acoustic local area network. Each array is able to simultaneously 
transmit different messages that will focus at the destination receiver node. Although the 
first protocols for acoustic local area networks have been proposed by Brandy and 
Catipovic (1994) and Talavage et al. (1994), the design of an acoustic local area network 
protocol that accommodates the learning step of the PCA’s is still yet to be developed. 

In order to achieve high-speed data transmission, the use of TRA bandwidth-efficient | 
phase-coherent communications requires further research. Due to temporal variability of 
channel conditions, the first step towards the feasibility of TRA phase-coherent 
communications is the evaluation of the phase-stability (phase fluctuations) at the focus. 
Once accomplished, the receiver structure can be determined and an adaptive equalizer 
can be designed to accommodate temporal variability of the channel conditions. 

At this point it seems clear that if the PCA learns every symbol in the code, the 
symbols will arrive with the correct phase information at the receiver. This has the 
downside that it will take much more time for the PCA to learn the environment and it 
requires enough storage for all of that information. It would be an advantage if the PCA 
only needed to learn the reference symbol and could then process the collected data in 


order to generate at the receiver location all the different symbols 1n the code. 
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